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1.0 INTRODUCTION 


The purpose of this research project is to apply control theory for dis- 
tributed systems to large flexible space structures, and to combine distri- 
buted system theory with finite dimensional control theory to produce 
comprehensive control system design methods for large space structures. The 
basic idea is to establish the existence and properties of an ideal infinite 
dimensional compensator for the distributed model of the structure and then 
approximate the ideal compensator with realizable finite dimensional compensa- 
tors. For an analysis of the ideal compensator and the resulting performance 
of the closed-loop system, distributed system theory, which is based on infin- 
ite dimensional analysis, is essential. Also, infinite dimensional analysis 
is necessary for studying convergence of approximations and selecting 
appropriate approximation schemes. For computing the approximating finite 
dimensional compensators and determining efficient implementation, finite 
dimensional control theory is used. 

To demonstrate the applicability of the methods to real space struc- 
tures, the research has focused on a model of a space antenna which consists 
of a rigid hub, flexible ribs and a mesh reflecting surface. This model was 
designed to include the most important characteristics of large flexible space 
structures. In particular, the structure to be controlled has different types 
of distributed components, as well as rigid-body modes, and, because of the 
mesh surface, there are many closely packed structural frequencies. This pro- 
ject appears to be the first application of distributed system control theory 
to a structure of such complexity. 
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It should be emphasized that, although the term "distributed system" 
usually implies a system involving partial differential equations (as with the 
model of the flexible components of the antenna here), the methods of this 
research dg not require the solution of partial differential equations. 

Rather, according to the approximation results, only finite dimensional Ric- 
cati and Liapunov matrix equations need be solved. 

Chapter 2 describes the space antenna model used in this project and 
discusses the finite element approximation of the distributed model. A Galer- 
kin component mode approach is used to develop a finite dimensional approxima- 
tion of the structure. This approach was found to converge rapidly to a model 
of suitable accuracy. 

The basic control problem is to design an optimal or near-optimal com- 
pensator to suppress the linear vibrations and rigid-body displacements of the 
structure, which is modeled as a distributed system because of the flexible 
components. As is common in finite dimensional optimal control of linear sys- 
tems, a quadratic performance index is defined to penalize both the distur- 
bance in the structure and the control effort used to eliminate the distur- 
bance. Determining the optimal compensator for this problem and analyzing the 
response of the system require linear-quadratic-gaussian (LQG) optimal control 
theory for distributed systems [BA1,CU1,GI1,2] . The application of this 
infinite dimensional LQG control theory to flexible structure is discussed in 
Chapter 3. 

Distributed system theory yields an infinite dimensional compensator as 
the optimal, or ideal solution to the control problem. As discussed in 
Chapter 3 this compensator can be represented in terms of functional gpins for 
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the feedback control and state estimator, and also in terms of an irrational 
transfer function. Both the functional gains and the Bode plots of the 
transfer function provide useful graphical representations of“the ideal com- 
pensator. 

Approximation theory is a central topic of this research. Because the 
ideal compensator is infinite dimensional, it can be neither written down in 
closed form nor built and implemented exactly. Thus, a finite dimensional 
approximation of the ideal compensator is required. In this research, there 
are two stages in the approximation of the ideal compensator. First, increas- 
ingly large-order finite dimensional LQG problems are solved until the result- 
ing compensators converge to a finite dimensional compensator that is essen- 
tially equal to the ideal compensator. Next, the order of this large compen- 
sator is reduced as much as possible using balanced realizations. It is 
important to note that, in both these stages of approximation, the convergence 
and truncation criteria are tied to quantitative measures of control system 
performance. 

The approximation theory for constructing the large-order approximation 
to the ideal compensator is discussed in Ref. [HR1] . At the heart of the 
infinite dimensional LQG theory, are the infinite dimensional Riccati operator 
equations whose solutions determine the optimal control and estimator gains 
for the infinite dimensional compensator. The scheme for computing finite 
dimensional approximations to the ideal compensator is based on the approxima- 
tion of the infinite dimensional Riccati equations by finite dimensional Ric- 
cati matrix equations of the type found in finite dimensional control theory. 
These finite dimensional equations are given in [HR1] . Also to be found in 
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[HR1] are typical functional control gains and the Bode plot for the ideal 
compensator for the antenna. 

While infinite dimensional, or distributed system control theory is the 
main basis for this research, finite dimensional control theory plays an 
important role. Indeed, the true richness of the project stems from the com- 
bination of a primary infinite dimensional perspective of control theory with 
finite dimensional control theory and practice. Thus, the project really 
represents an application of a finite dimensional methods in a infinite dimen- 
sional setting. 

After the converged approximation to the ideal compensator has been 
obtained, it must be verified that the lower order approximation to the ideal 
compensator will still produce a satisfactory closed-loop system. Also, 
robustness with respect to errors in plant parameters must be studied. 
Robustness issues were a major focus of the work performed this year. Two 
basic approaches for robustness enhancement were investigated: Loop Transfer 

Recovery and Sensitivity Optimization. A third approach synthesized from ele- 
ments of these two basic approaches is currently under development. The goal 
of loop recovery is to produce a loop transfer matrix with the same frequency 
response as full state feedback which is known to have very desirable robust- 
ness properties. Sensitivity optimization uses numerical optimization tech- 
niques to select regulator and/or estimator gains which reduce sensitivity of 
the real parts of the closed loop poles to parameter errors. These methods 
are presented in Sections 4.2 and 4.3 along with numerical examples which 
illustrate their performance. 
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In general it was found that the methods worked well on relatively sim- 
ple problems which were used to test and develop the software. However, 
results were mixed when the methods vere applied to the antenna model. The 
light damping and complexity of the antenna apparently make this a challenging 
system for robust controller design. The performance requirements were pur- 
posely chosen to ensure that several of the flexible modes were within the 
controller bandwidth, and this contributed to the difficulty of the problem. 
The form of the weighting matrix Q in the performance term of the cost func- 
tion was found to have a strong effect on robustness. For example, Q's which 
did not penalize velocity terms were found to be less robust than those which 
did. Also, Q's in which the penalty on rigid body motion was much greater 
than the penalty on flexible modes were less robust than Q's in which the 
rigid body mode was not emphasized. These results are reported in Section 
4.2. 


Sensitivity optimization, a second method for designing robust compensa- 
tors, is discussed in Section 4.3. This method uses nonlinear programming to 
find control and estimator gains that minimize the sensitivity of the closed- 
loop response with respect to uncertain plant parameters. Of the various 
measures of closed-loop sensitivity that were tried, the sensitivity of the 
real parts of the closed-loop eigenvalues proved to have the best correlation 
to robustness. Simulation results are given. As with the loop-recovery 
methods, the form of state weighting in the control problem affects the 
robustness significantly. 

Chapter 5 discusses control driven finite element approximation of flex- 
ible structures. This means choosing the finite element scheme that is most 
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efficient for approximating the solution to the optimal control problem for 
the distributed model of the structure* Section 5.1 compares three sets of 
finite element basis vectors for computing functional control gains and demon- 
strates that some finite element schemes are much better than others. Section 
5.2 discusses the possibility of constructing a finite element scheme to 
approximate the infinite dimensional Hamiltonian system directly, instead of 
indirectly as in Chapter 3. The main application of this approach appears to 
be for structures with history -dependent damping. 
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2 


STRUCTURE MODELS AND FINITE ELEMENT APPROXIMATION 


2.1 ANTENNA MODEL 

The space antenna model that has served as the primary example of this 
research is shown in Figure 2.1.1. The flat antenna consists of a rigid hub, 
eight ribs and a mesh reflecting surface. The ribs are modeled as beams can- 
tilevered to the hub, and the mesh is modeled as circular sectors of membrane 
tied to the ribs and hub. The center of the hub is fixed but the hub can 
rotate out of plane, so that there are two rigid-body modes. 

This model is based on the Lockheed wrap- rib antenna [L01] . Since this 
is the first complex structure to which the methods of this research have been 
applied, the model has been simplified to a flat antenna with eight instead of 
the actual 48 ribs, to allow concentration on fundamental issues while main- 
taining a complex structure with different types of components. Otherwise, 
the parameters given in Figure 2.1.1 are based on the 48-rib antenna [ELI]. 

The two actuators apply torques to the hub, and the sensors measure the 
rotation of the hub and the displacement of the tip of each rib. The compen- 
sator is designed to control the out-of- plane motion of the antenna. 

For small elastic displacements, the symmetry of the antenna further 
reduces the complexity of the compensator because the motion of the antenna 
decouples into two sets of orthogonal modes, with each set controlled indepen- 
dently by one actuator. Each of these sets consists of modes that are asym- 
metric about one torque axis. Although the actuators can control only the 
controllable modes, these are the only modes excited by rotating the antenna 
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with the torques on the hub. 


The equations of motion for rigid-body rotation and small elastic vibra- 
tion have the form discussed In Chapter 3. See (3.1.1) - (3.1.10). Here, we 
will discuss the stiffness and damping operators, which are the primary opera- 
tors In the equations of motion. The stiffness operator A q can be separated 
into rib and mesh components as 

4 o-V + W (2.1.1) 

where 


A Qr - El d 4 (.)/dr 4 . 


( 2 . 1 . 2 ) 


A om - Tild 2 <.)/dr 2 + (l/r)d(.)/drl - T, d 2 <.)/dd 2 

(2.1.3) 

and Tj and T 2 are the tensions in the radial and transverse directions, 
respectively. The term A Qr operates on the elastic deformation of each rib 
an< * ^om °P erates on the elastic deformation of the mesh. 

We assume that both the ribs and the mesh have viscoelastic damping, so 
that the damping operator is 


D ° ■ * + * (2.1.4) 

where Q p and Cj, are positive scalars. If c r and Cj, are not equal, this struc- 
tural damping couples the modes — but not modes that are controllable by one 
actuator with modes that are uncontrollable by the same actuator. 


8 



MODEL DATA (based on 48 rib Lockheed wrap rib antenna): 

Hub Radius: 46 in. Hub Weight: lOOO lbs. 

Hub Inertias: I MM = I vv — 342 lb— in— sec* 

I«. = 684 lb-in-sec* 

Rib Length: 86 -ft. Rib Weight: 115 lb. per rib 

Rib sti-f-fness: El *» 4.05xl0' r lb— in* 

Mesh Weight: 6.94 lb per sector 

Figure 2.1.1 Q 


2.2 FINITE ELEMENT APPROXIMATION OF THE ANTENNA 


The approximation of the antenna is a component modal synthesis for 
which the component modes of the ribs and mesh were determined by consistent- 
mass finite element approximations. The ribs were approximated by a standard 
finite element scheme using Hermite cubic splines, and the mesh sectors were 
approximated by a Galerkin scheme discussed in Section 2.3. 

In an earlier project [HR1] , the approximation scheme used in this pro- 
ject was compared to a NASTRAN finite element approximation of the antenna. 
While both schemes converged to the same mode shapes and frequencies as the 
order of approximation increased, the component modal synthesis converged fas- 
ter and was more convenient and cheaper to use. An important advantage of the 
conponent modal approximation for this project is the fact that the program 
was written to make the order of approximation easy to vary. In design of 
controllers for complex structures, it is impossible to determine in advance 
what order of approximation will be necessary for designing an optimal, or 
even acceptable compensator. Thus, in a project like this, it is essential to 
have an approximation scheme of easily variable order. 


2.3 GALERKIN APPROXIMATION OF SECTOR 


A Galerkin approximation of the antenna was used to compute the optimal 
control law for the composite sector. The basis vectors have the form 



(2.3.1) 


where (0 n< e£) are the natural modes of the hub and ribs without the mesh, and 
<P^ are the natural modes of the mesh when the beams and hub are constrained 
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not to move. This Galerkin approximation Is then a component modal synthesis 
of the composite sector. 


The mesh modes have the form 

v“(d» r) * sin8nd f n# k(r), n - 1, 2, .... 


k * 1, 2, .... 

where the functions ^(r) are the eigenfunctions of the problem 

T 1 [8(64n 2 /r 2 )f - (f” + f’/r)] = “J,k f 


(2.3.2) 


(2.3.3) 


f(r o ) - r (r 0 + l r ) - 0 


(2.3.4) 


Here, the eigenvalues « n ^ are the squares of the natural frequencies of the 
mesh, i s the radial tension, Tj/Tj ■ 0 = 72 and & r is the rib length. This 

value of 0 was selected so that the first mesh frequency is approximately half 
way between the first two controllable beam-hub frequencies. 


A standard Ritz-Galerkin approximation with linear splines as basis vec- 
tors was used to solve (2.3.3) and (2.3.4) for the mesh modes. As many as 43 
splines were used, but convergence of the first 10 frequencies and mode shapes 
for n = 1 was obtained with 30 splines. 
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3.0 STATEMENT AND SOLUTION OF CONTROL AND ESTIMATION PROBLEMS 


3.1 ABSTRACT CONTROL SYSTEM 


The space antenna model used here falls Into the large class of flexible 
structures represented by the abstract control system 



V 


At 

o 


B 0 u + (/p 


(3.1.1) 


y 


C o x 


+ </« 


(3.1.2) 


where the generalized displacement x(t) is in a real Hilbert space H, the con- 
trol u(t)« R®, the measurement y(t)« R p . The disturbance (4 and the noise U 0 
are zero-mean Gaussian white noise processes, as in the finite dimensional LOG 
problem. The mass operator M q i S bounded, selfadjoint and coercive on H, B Q 
is bounded from R® to H and C Q j. s bounded from H to R p . The unbounded stiff- 
ness operator A q i S selfadjoint and bounded below, with compact resolvent, and 
the damping operator D q i a symmetric, nonnegative and A Q - bounded. A similar 
abstract description of the LQG control problem has been used by Balas in 
[BA2] and [BA3] and Gibson in [Gil]. Also, see Balakrishnan [BA 1 ] . 


As usual, by natural modes we mean the eigenvectors Xj 0 f the problem 


M_ 


{ j * A o x j 


1 . 2 , 


(3.1.3) 


where the are the natural frequencies of the structure. In the antenna 
problem, as in most structure control problems, the only nonpositive eigen- 
values of A q are zero eigenvalues corresponding to rigid body modes. When 
these rigid body modes are controllable, they are contained in the positive 
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eigenspace of B q b*, so that the operator 


A o “ A 0 + B q B* 2 m 2 0, 


for some positive m; i.e. , A q i S coercive. 


Next, we define the strain-energy space V * D(A^ 2 


<V l' y 2>V " <A o 2v 1. ^ 2y 2 > H’ 


(3.1.4) 

) with inner product 

(3.1.5) 


This is a natural space for the generalized displacement because usu- 
ally, as in our example, <x,x> y is just the sum of twice the elastic strain 
energy and the squares of the rigid body displacements. 


To write (3.1.1) in first order form, we define the energy space E - 
VxH M where H m is H with the equivalent kinetic-energy inner product 


* <Mhj,h2>n. 

• 2 

Thus || xll M is twice the kinetic energy in the structure. 
(3.1.1) as 


(3.1.6) 

Now we can write 


z ■ Az + Bu + (/,, 

(3.1.7) 


y 


Cz 


+ </ 
w o* 


(3.1.8) 


where 


z ■ (?) 
x 


0 


0 


B - 

M _1 B 
n o o 


(3.1.9) 


and 
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A 


D(A 0 ) x D(A 0 )CD(A) 


-M^A ~X~ n \ 

O O 0 0 


(3.1.10) 


As In [Gil, Sec. 2], D(A) Is chosen uniquely so that A generates a C Q - 
semlgroup T(t) on E. 


3 .2 INFINITE DIMENSIONAL LQG PROBLEM 


As in the finite dimensional case, a separation principle [BA1, CU1] 
allows us to design the optimal compensator by solving the deterministic 
optimal regulator problem and the stochastic state-estimation problem 
separately. The infinite dimensional optimal control problem is to choose u 
to minimize 


J * y (<Qz, z> E + <Ru,u>jjm)dt, 

0 (3.2.1) 

where z is the solution to (3.1.7), Q = Q* is a nonnegative bounded linear 

operator on E, and H * R*>0. 

To simplify certain technicalities about existence of solutions to the 

LQG problem, we will assume that + b„B* is coercive, which means that the 

o o o 

flexible components of the structure have coercive damping, and that Q is 
coercive, which guarantees that the optimal closed-loop system is uniformly 
exponentially stable. 

For this deterministic problem, the optimal control has the feedback 

form 


u » — Kz 



(3.2.2) 


where 


K - -R _1 Bll 


(3.2.3) 


and n is the unique nonnegative selfadjoint element of L(E) which satisfies 
the Infinite dimensional Riccati equation 


A^I + DA - DBR _1 B* n + Q « 0 


(3.2.4) 


See [BA1. CU1, Gil, G12]. 


The minimum variance estimator, or infinite dimensional Kalman filter. 


is [BA1. C02] 


z * Az + Bu + G(y-C&) . 


(3.2.5) 


where 


A i 

G = IKR -1 


and II satisfies the Riccati equation 


A . A *A . A a 

JU + DA - IK R -1 C D + Q 


(3.2.6) 


(3.2.7) 


A A A* 

The pxp matrix R and the bounded nonnegative operator Q =* Q are the 
covariance operators for (/ Q and (/*, respectively. We assume that any undamped 
modes are observable, so that the closed-loop estimator is uniformly exponen- 


tially stable. 
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We do not require either Q or Q to be trace class, as Is usually neces- 
sary for the Infinite dimensional LQG problem [BA1, CU2] . This Is not neces- 
sary for our definition of the Infinite dimensional compensator, and It allows 
more freedeom In designing the compensator to produce desired closed-loop 
transient response. See also [GI2] . 


The optimal compensator consists of the Infinite dimensional estimator 
(3.2.5) and the control law 



where K Is given by (3.2.3) with (3.2.4), and (3.2.5) becomes 

A a A 
z m A c z + Gy, 


where 


(3.2.8) 


(3.2.9) 


A 


A - 


- GC 


(3.2.10) 


Note that the optimal compensator has the Irrational transfer function 


KlsI - A] -1 G 
c 


(3.2.11) 


Even though this transfer function would require an infinite dimensional 
state-space realization, the transfer function itself is only a finite dimen- 
sional matrix function of the complex variable s. For control of the antenna 
quadrant, this Is just a 1 x 3 row matrix representing the three control chan- 
nels from the displacement sensors at the rib tips to the torque actuator on 
the hub. 
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3.3 FUNCTIONAL GAINS 


When there is a single actuator (nr>l). the operator B Q i s actually an 
element of H and BsG. Hence 


or 


Kz - R _1 b1Iz « flIBR _1 ,z> E , 


(3.3.1) 


where 


Kz * <f,x> v + <g,i> M 


(3.3.2) 


(£) -DBR" 1 

o 

with feV and geHjj - H. We call f and g functional control gains . 


(3.3.3) 


Similarly, for a single sensor, the measurement y is a scalar and C*eE. 


Then 


G 


A * i 

ncR 1 


A 



(3.3.4) 


A A 

where feV and geH are functional estimator gains . 

For the multi-input-multi-output case, there is a pair of functional 
gains for each actuator and for each sensor. 


3 .4 APPROXIMATION 

To obtain a sequence of approximating finite dimensional LQG problems, 
we use a Ritz-Galerkin approximation scheme for (3.1.1). We assume a sequence 
of linearly independent basis vectors e,, which are complete in V. The n th 
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approximate solution to (3.1.1) Is 


V fc > “ ^*i (t>e i' — 

1-1 (3.4.1) 

which is in « span (e 1# ,e Q ). * e need both V n and H Mn , which are 

the same set but have the V and inner products, respectively. 

Initially, let us oonsider the case (/j - 0 in (3.1.1). Then the coeffi- 
cients a n (t) satisfy 


M V” + D I 11 + A a n ■ B 
on a on a on a 


on 


u. 


where a n (t) is the n-vector containing a n (t). i ■ 1, .... n. 
ing and stiffness matrices are given respectively by 


(3.4.2) 
The mass, damp- 


en, ij * <e i ,e j > M " <M o e i' e j > H' 


(3.4.3) 


D on, ij = < D o e i' e j > H' 


(3.4.4) 


A on,iJ “ < A o e i' e j>H < 3 . 4.5 

In general. (3.4.4) and (3.4.5) are valid only if the basis vectors are 
in D < A 0 > * Otherwise, we use 

A on.ij “ <e l’ e J > V “ <B o®i' B o e j > H (345 

(Recall (3.1.4) and (3.1.3)) and a similar expression for D qq (i n our model, 
is essentially a scalar times A Q .) Also, 
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B 


on, ij 


<e i ,B 0 ej>, i - 1, 


J “ 1» 


(3.4.6) 


The convergence of such approximations and the corresponding finite 
dimensional optimal control problems is discussed in [3.5] for e ie D(A Q ), which 
includes the cases where the basis vectors are either natural modes of the 
structure or component modes. Here, we only outline the formulation of the 
sequence of n-order LQG problems. The most efficient way to do this is to 
note that, with (3.4.1), (3.4.2) is equivalent to the following differential 

equation on - span (ej e n ) : 

• • • 

M on x n + D on x n + A on x n “ B on u * (3 4 7) 

where M Qn , D on , A on , and B on are the operators on determined by the 
matrices in (3.4.3) - (3.4.5) and the identification (3.4.1). Of course, we 
can write (3.4.7) as 


n 


Vn 


B n u 


with z. 


(x 


n* 


x n> * 


(3.4.8) 


For (3.4.7) and (3.4.8), the optimal regulator problem leads to the fin- 
ite dimensional Biccati equation 


A r^n + ^ A n “ ^n B n R B 5^n + Q n * 4 9) 

where the operator Q n ± a defined as follows. Let E n - V n x and denote by 
P Q the projection of E onto E Q . Then 

Q n “ WilE • 

" (3.4.10) 

With this Q n and our preceding hypotheses about damping and the completeness 
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of the e^'s i n v, the 11^ of (3.4.9) is guaranteed to converge strongly to the 
solution II of the Infinite dimensional Riccati equation (3.2.4). 


The functional control gains in (3.3.3) are approximated by 


<«“> -W 1 - 

(3.4.11) 

The strong convergence of 13^ and B n implies that f n converges in V to f and g n 
converges in H to g. 

The functional control gains f and g n are associated with the n th order 
control law 


u - - K 


n z n " “ ^n' x n*V " ^®n ,x n^H* 


(3.4.12) 


with 


S n - »' lB X 


(3.4.13) 


To approximate the Infinite dimensional compensator, we construct a fin- 
ite dimensional state estimator 


where 


*n - A n z n + V + G n ( y " c n z n > • 




(3.4.14) 


G n mU n c n R 
A 

and ^n satisfies the Riccati equation 

A A | A ^ A ^ 

A A 'Vo ’ C fln ♦ »n - «• 


(3.4.13) 


(3.4.16) 
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A 

The operator Q q is given by 


A 


A 

P n GP n 


(3.4.17) 


for Q the covariance of the process noise <J X in (3.1.7). As before, R is the 
covariance of the measurement noise U Q . The operator P n is the projection 
onto E n . y n x Hjfo. 


He now have the components of the n th approximation to the infinite 
dimensional compensator. The n^ finite dimensional compensator oonsists of 
(3.4.13) and the oontrol law 


This compensator has the rational transfer function 

* n I=I - A*,]' 1 C n 

where 


(3.4.18) 


(3.4.19) 


cn 


A n - B n K n " 


G n c n* 


(3.4.20) 


For each sep(A c ), this transfer function (see Section 9 of [GI2]) the value of 
the transfer function in (3.2.11) as n Increases. 


3.5 APPLICATION TO SPACE ANTENNA 

The generalized displacement vector x for this system has components 
representing each of the two out-of-plane rigid body angles, the out-of-plane 
elastic deformation of each rib and of each mesh. The basic space H is then 
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H - R 2 x L 2 (0,1) x L 2 (0) x 


X L 2 (0,1) * L 2 (0) 


(3.5.1) 


for each rib and mesh 

where 1 is the length of the ribs and Q is the area of each mesh sector. The 
strain-energy space V is like H with L 2 (o,l) x L 2 (Q) replaced by H^O.l) x 
H*(fl) . 

In the numerical studies in [HR1] , we used the following data. For the 
performance index in (3.2.1), the state weighting operator Q and control 
weighting matrix R are taken as 


Q ■ I and R ■ .1x1. 

(3.5.2) 

Hence <Qz,z>g is twice the total energy in the antenna, plus the sum of the 
squares of the rigid body angles. 


The disturbance (J^ in (3.1.7) and the measurement noise (/ 0 in (3.1.8) 
are assumed to have the respective covariances 


A 

Q - 


[$ 


°J 


A 

said R * .01 x I. 


The viscoelastic damping coefficients in (2.1.2) are 


(3.5.3) 


'j, * .001 and Cj, = .003, 


The functional gains have the form 


f " («i» o 2 ,0a» d?* •••» d?) 

functions for each rib and mesh 


(3.5.4) 


(3.5.5) 
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4.0 ROBUSTNESS 


4.1 OVERVIEW 

Robustness refers to the ability of a control system to perform satisfacto- 
rily even when the model used for design is an imperfect representation of the 
actual system. Modeling difficulties are a major concern in designing control 
systems for LSS. Much of the work on functional gains , described in the 
Ref. [HR1] , is aimed at one particular type of modeling problem: model trun- 

cation. Parameter errors are another important source of modeling error. 
During the current year of research, we have worked on the problem of develop- 
ing LSS controller designs which are relatively insensitive to parameter 
errors, but still maintain a satisfactory level of system performance. For 
antenna structures of the type considered, this is a difficult problem which 
does not easily yield to standard techniques. Because of its importance, how- 
ever, it is a problem which cannot be neglected. 

Two different techniques for robustness enhancement have been explored dur- 
ing 1985: Loop Transfer Recovery (LTR) and Sensitivity Optimization. Sec- 

tions 4.2 and 4.3, which follow, describe the basic features of these 
approaches, and also give examples of their performance. Section 4.4 briefly 
discusses a controller design approach based on elements of both Loop Transfer 
Recovery and Sensitivity Optimization. At the present time this synthesized 
approach is still under development. 

> 

As a general rule, an improvement in robustness can only be achieved at the 
cost of a reduction in some measure of performance. In order to achieve a 
proper balance between performance and robustness, some constraints must be 
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imposed. In the loop recovery approach, the open- loop and closed- loop 
frequency responses are constrained to approximate those of the full state 
feedback LQR design. Thus, the performance of the system is first determined 
in the context of the LQR design, and then one tries to obtain an estimated 
state feedback implementation which has similar performance and' also adequate 
robustness. The sensitivity optimization approach also begins with the solu- 
tion of an LQR problem and designs an estimated state feedback implementation. 
In this case, however, there is no constraint to retain the LQR frequency res- 
ponse. Instead, the goal is to minimize the sensitivity of the closed-loop 
eigenvalues to parameter errors. A constraint must be imposed to keep the 
closed loop regulator and estimator eigenvalues sufficiently far in the left- 
half -plane, otherwise the optimization routine would generate a very low per- 
formance control system. Designs based on sensitivity optimization tend to be 
more robust than LTR designs, but their loop gains are also lower. 

Both of these approaches rely heavily on the statement of the LQR problem 
which serves as a starting point for the designs. The state weighting opera- 
tor in the performance index for the LQR problem discussed in Chapter 3 is Q. 
Generally, scaling up Q while holding the control weighting matrix R fixed 
leads to larger control gains and a less robust compensator. The numerical 
examples presented in Sections 4.2 and 4.3 indicate that in addition to the 
magnitude of Q, the relative penalization that Q specifies among the different 
components of the state vector can affect robustness dramatically. For exam- 
ple, Q's which penalize only the RMS 1 surface error of an antenna model were 

1 Throughout this report, the term RMS surface error refers to the root mean 
square displacement of the antenna surface relative to its nominal position 
in inertial space. Thus RMS surface error includes the effects of rigid 
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found to lead to designs which were not very robust. However, Q's which 
penalize the energy or the RMS surface error plus the rate of change of RMS 
error, were found to be quite robust. While the interpretation of these 
results is not entirely clear at this point, some observations can be made. 
Penalizing RMS surface error alone leads to controllers which try to reduce 
rigid body displacement much faster than elastic deformations. This means 
that the closed loop eigenvalues for rigid body modes are substantially to the 
left of those for the lowest flexible modes. In this case, large control tor- 
ques are produced to control rigid body modes, and the system is not robust 
because these torques strongly excite the flexible modes, which may be inaccu- 
rately known. The other two Q's mentioned above penalize rigid body and flex- 
ible modes more equally. The real parts of the closed-loop eigenvalues, and 
hence the decay rates, for rigid body modes and lowest elastic modes are 
approximately equal and the compensator is more robust. In our continuing 
research effort, we plan to pursue further the relationship between the selec- 
tion of Q and robustness. 


body rotation as well as elastic deformations of the antenna 
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4.2 LOOP TRANSFER RECOVERY TECHNIQUES 


It is well known that linear quadratic regulator designs have very attrac- 
tive robustness properties, e.g., one half to infinite gain margin and 60° 
phase margin. This property is based on the use of full-state feedback. In 
practice, full-state feedback is virtually never available (certainly not for 
LSS) and the state must be estimated. The inclusion of an estimator fre- 
quently results in loss of the robustness promised by the regulator design. 
Loop recovery techniques have been proposed as a way to recover a measure of 
this robustness [D01,2]. The Loop Transfer Recovery (LTR) approach was origi- 
nally suggested by Kwakernaak [KW1] and later extended by Doyle and Stein 
[D01,2]. Kwakernaak first derived a method to asymptotically recover the loop 
shape of a given Kalman-Bucy Filter (KBF) at the plant output, by allowing the 
control weighting in a specific LQR problem to approach zero. This offers the 
advantage of recovering the minimal sensitivity properties of the Kalman Fil- 
ter. Doyle and Stein pointed out that the dual to Kwakernaak 1 s approach, 
allowing the noise covariance in a KBF problem to approach zero, asymptoti- 
cally recovers the LQR loop shape at the plant input. This offers the advan- 
tage of recovering the gain and phase margins of the Full-State Feedback Quad- 
ratic Regulator. Both procedures will be referred to as Loop Transfer 
Recovery (LTR), where Kwakernaak 's approach [Ktfl] recovers a loop shape at the 
plant output, and is only valid when the plant has at least as many inputs as 
outputs, while the approach of Doyle and Stein [D01,2] recovers a loop shape 
at the plant input, and is only valid when the plant has at least as many out- 
puts as inputs . Both approaches require that the plant be non-minimum phase 
and it is only in the case of square plants that the designer has a choice of 
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recovering the loop shape at either the input or output 1 . A recent paper 
which discusses more fully the various aspects of examining the feedback loop 
at various points is Ref. [FR1] 

The present work examines the application of the LTR approach to designing 
robust control systems for lightly damped, flexible structures with an unequal 
number of inputs and outputs. It is relatively common for a flexible struc- 
ture to have more available measurements than available control inputs. 

Hence, the case of an excess of outputs over inputs will be emphasized. In 
particular, the special case of 1-input and m-outputs, which, by reducing the 
loop gain at the plant input to a scalar function, both simplifies the analy- 
sis and clarifies conceptual ideas, will be used extensively. This is moti- 
vated by the study of the wrap-rib antenna. 

The first implication of the excess of outputs over inputs is that the loop 
shape can be recovered only at the plant input, while a second implication is 
that a number of extra degrees of freedom, required to produce a square sys- 
tem, are available. The implication of these degrees of freedom will be care- 
fully considered and two methods for taking advantage of them will be offered. 
In particular, an algebraic method that will provide Loop Transfer Recovery 
with arbitrary compensator pole placement and a reduced order compensator, 
will be presented. An optimization method which chooses the extra degrees of 
freedom so as to minimize the sensitivity of the closed-loop pole locations 
with respect to plant frequency errors is discussed in Section 4.4. 


1 Margins at both the input and the output can be guaranteed simultaneously 
only when the inequalities 7[I+K(s)G(s) ]kl and 7T[I+G(s)K(s) ]£1 both hold 
[LEI]. 
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It should be noted that at least one paper on an application of LTR to a 
flexible space antenna has appeared in the literature [SU1] . In this work the 
design model was reduced to the three rigid body modes of the antenna, while 
all the flexible modes were treated as uncertainties. Unstructured uncertain- 
ties are a good representation of the error due to neglected dynamics so this 
is a valid application of LTR, but it will clearly result in a very low per- 
formance control law 2 . In fact, the resulting control law indicates an open- 

-3 

loop bandwidth below 10 rad. /sec., where the first flexible mode is not 
encountered until approximately .75 rad. /sec. The application studied in the 
current work involves finding high performance control laws for flexible 
structures with uncertain frequencies. This means that the control action 
will actively control the uncertain frequencies, or equivalently, open-loop 
bandwidths will fall beyond the uncertain frequencies. In this case the unc- 
ertainties are best described in terms of parameter errors rather than unmo- 
deled (i.e. truncated ) dynamics. As illustrated in Section 4.2.1, unstruc- 
tured uncertainties are not a good representation of parameter errors in the 
modeled modes of a lightly damped flexible system and the difficulty in 
achieving a robust design is greatly increased. 

The organization of the next few Sections is as follows. First the ques- 
tion of robustness is discussed, paying particular attention to the modern 
robustness Theorems of the past 10 years. The most important results relevant 
to this study are the unstructured representations of plant uncertainty pre- 
sented by Doyle and Stein [D01] and by Lehtomaki [LEI, 2, 3, 4]. The unstruc- 

2 The reason why this implies a very low performance control law is discussed 
in Section 4.2.3. 
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tured uncertainty indicates that robustness can be guaranteed by making the 
minimum singular value of either the return difference or inverse return dif- 
ference transfer function matrices large enough. This leads to the idea of 
defining robustness in terms of loop shape, and provides the motivation for 
loop shaping as a control design method, and LTR as a particular approach to 
achieving some desired loop shape. Using a simple one mode example, it is 
shown that the unstructured uncertainty is far too conservative to treat 
poorly modeled frequencies. It is concluded that in the case of a strict par- 
ameter uncertainty, such as an uncertain frequency, robustness should be 
checked by varying the parameter within the expected range, while examining 
the closed loop eigenvalues. A secondary, but crucial conclusion is that 
robustness is not necessarily determined by loop shape, except in the special 
case where the system component uncertainties are accurately described in 
terms of unstructured uncertainties . 

Next a description of the LTR method is given, along with a discussion of 
some of its properties. Again, particular attention is paid to the implica- 
tions of some of these properties for lightly damped, flexible structures. It 
is concluded that the LTR method does not always guarantee a robust system, 
even when the corresponding full-state feedback design is extremely robust, 
and should be used carefully in designing compensators for lightly damped, 
flexible systems. This is a result of the previously noted conclusion, that 
robustness is not necessarily a function of loop shape alone. The reasons for 
lack of robustness in LTR designs is examined and methods for using the extra 
degrees of freedom available in the non-square problem, to find a "better" LTR 
design for a given LQR loop shape are suggested. 
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After presenting the design methods, they are applied to a wrap-rib antenna 
model, in an attempt to achieve robustness in the face of frequency uncer- 
tainty, while maintaining a given loop shape. The antenna model is chosen 
because it is much more sensitive to frequency errors than other models we 
have tried. The results are compared with traditional LQG designs, tradi- 
tional LTR designs and with a direct parameter optimization approach. In the 
conclusions, the implications of this comparison, and a general approach to 
using LTR to design robust control systems for non-square, lightly damped, 
flexible structures are discussed. 

Background information may be found in the appendices, including a review 
of singular values, proofs of Theorems, and a discussion of numerical consid- 
erations in the application of the two LTR methods. 


4.2.1 ROBUSTNESS MEASURES FOR MIMO FLEXIBLE SYSTEMS 

One of the major reasons for the application of feedback control, is to 
minimize the effect of variations in plant dynamics (or equivalently, plant 
frequency response) on the system performance. Robustness is the study of 
exactly how large of a variation in plant dynamics can be tolerated, before a 
given feedback system goes unstable. 

The basic work in single-input-single-output (SISO) feedback systems was 
done more than forty years ago by Nyquist, Bode [B01] and their colleagues. 

The most important robustness result is the Nyquist Stability Criterion, and 
the related concepts of gain margin and phase margin. These margins specify, 
in an exact sense, how much the gain and phase of the plant frequency response 
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can vary before the closed- loop system will go unstable. In particular the 
gain margin measures the amount by which the open- loop gain can be increased 
without causing instability, while the phase margin measures the amount by 
which the open- loop phase lag can be increased without causing instability. 
These essentially measure the closeness of the Nyquist plot to the -1 point, 
therefore generalizing the Nyquist Stability Criterion from a measure of abso- 
lute stability to one of relative stability. 

The present work focuses on developing robust control methods for large, 
flexible, lightly damped structures. There are two aspects of these systems 
that cause some difficulty in analyzing robustness properties. The first is 
their lightly damped, flexible nature, which implies a highly oscillatory fre- 
quency response, making concepts such as gain and phase margin much more dif- 
ficult to interpret. This difficulty, however, can still be handled within 
the context of classical control methods, though with some care. The second 
difficulty proves to be more fundamental, and stems from the existence of a 
number of inputs and outputs, all connected together. Since there is no lon- 
ger one loop gain to be analyzed, classical methods cannot be applied 
directly. Furthermore, it has been shown that analyzing each loop of a multi- 
ple loop system separately, does not always give results that are valid for 
the overall system [ DO 1 ] . Fortunately the robustness of multivariable systems 
has been studied extensively in the last ten years. It is the purpose of this 
Section, to give a quick overview of the state-of-the-art in robustness of 
multiple-input-multiple-output (MIMO) feedback control systems. 

An extremely general approach to MIMO robustness can be found in the work 
of Zames [ZA1,2] on cone bounded perturbations and Safonov [SA1,2] on even 
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more general perturbations. These results rely heavily on the mathematical 
field of functional analysis and are applicable to both non-linear and time- 
varying plants. The special case of linear, time -invariant plants, however, 
can be dealt with much more simply, and a number of more practical robustness 
Theorems have been developed. These Theorems can be divided into two groups, 
paralleling the two modem approaches to control theory mentioned in the last 
Section. The most popular approaches to describing the robustness of linear, 
time- invar iant MIMO feedback control systems deal with a transfer function 
description of the plant uncertainty. It's the transfer function approach 
which leads to loop shaping as a method for obtaining robust control, and it 
will be dealt with first. 


Transfer Function Approaches to Robustness 


Consider the simple feedback system illustrated in Fig. 4.1, where: 

R(s) - Command Signal 

K(s) - Compensator 

U(s) - Control signal to the plant 

G (s)- True, possibly unknown plant, described by a nominal 
plant, G(s), and some characterization of errors 
Y(s) - Output variables (Measurements) 

Dropping the dependence on s, some transfer functions describing the system 
can then be evaluated as follows: 


Y = [ (I+GK) -1 GK]R = [GK(I+GK) _1 ]R 
E = R - Y = (I+GK) _1 R 
U = [ (I+KG)” 1 K]R = [K(I+KG) _1 ]R 

Also define the following functions: 


( 4 . 2 . 1 ) 

( 4 . 2 . 2 ) 

( 4 . 2 . 3 ) 
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6K - Output Loop Transfer Function 

KG - Input Loop Transfer Function 

[I+GK] - Output Return Difference 

[I+KG] - Input Return Difference 

[I+(GK) *] - Output Inverse Return Difference 

[I+(KG) - Input Inverse Return Difference 

And make note of the following identities: 


[I+CGK)" 1 ]' 1 = GK[I+GK] _1 = [I+GK] _1 GK (4.2.4a) 
[I+CKG)" 1 ]” 1 = KG [I+KG] -1 = [I+KG] _1 KG (4.2.4b) 
[I+(GK) -1 ] _1 + [I+GK] -1 = I (4.2.5a) 
(I+(KG) -1 ] _1 + [I+KG]" 1 = I (4.2.5b) 


The first two identities show that the inverse return difference is simply the 
inverse of the closed loop transfer function. The second two identities show 
that the sum of the inverses of the return difference and inverse return dif- 
ference is equal to the identity. The implications of this fact, in terms of 
fundamental limitations on design, are discussed in Ref. [SA3] . 

Note that the output loop and input loop transfer functions and related 
return differences are not in general the same for MIMO systems, though they 
are for SISO systems. In fact, if the number of inputs is not equal to the 
number of outputs, they will not even have the same dimensions. This property 
affects the way in which an uncertainty is described. In particular the true 
plant G'(s) might be described in terms of a nominal plant G(s) in one the 
following two ways: 
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(4.2.6a) 


G'(s) = L o (s)G(s) 

-or-— G'(s) = G(s)L i (s) (4.2.6b) 

L q (s) will have the dimension of the plant output, and can be considered an 

uncertainty acting at the output, while L^(s) will have the dimension of the 

plant input, and can be considered an uncertainty acting at the input. L q (s) 

would provide the best error description for the output loop since 
G' (s)K(s)=L q (s)G(s)K(s) , while L^(s) would provide the best error description 

for the input loop, since K(s)G' (s)=K(s)G(s)L^(s) . The point at which the 

loop is "opened" to examine robustness will therefore depend on the way in 
which the error is described, and vice-versa. For simplicity, only errors 
acting at the plant input will be considered, but the results for errors act- 
ing at the output can be easily found by replacing K(s)G(s) by G(s)K(s) in the 
Theorems of the following Sections. 


Unstructured Uncertainty 

The most basic description of uncertainties in the form of (4.2.6a) or 
(4.2.6b) is the unstructured uncertainty introduced by Doyle and Stein [D01] 
and by Lehtomaki [LEI, 2, 3, 4] . An unstructured uncertainty is one that can be 
given no more structure than a bound on a suitable measure of its size. For 
matrices, singular values serve as such a measure. Singular values are dis- 
cussed in Appendix A. It is assumed that the reader is familiar with the con- 
cept of maximum and minimum singular values as well as some of their proper- 
ties. 
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Consider the following two 1 possibilities for L(s) = L^s): 

L(s) = [I+A^Cs)] (multiplicative uncertainty) (4.2.7a) 

or L(s) * (I+A^(s)] * (divided uncertainty) (4.2.7b) 

where in each case the perturbation can be bounded by frequency dependent 

functions 1 and l.: 2 
m d 

■o[A m (s)] £ l m (s) (4.2.8a) 

and '3'[A < j(s) ] £ l d (s) (4.2.8b) 

THe function l ffl (s) and l d (s) are real and positive on the Nyquist D-contour. 

The closed loop stability of the true system (KG'[I+KG'] *) can be deter- 
mined by the multivariable generalization of the Nyquist Criterion [R01], 
which requires that the det[I+KG'], evaluated on the standard Nyquist D-cont- 
our (denoted by seft^), encircle the origin in a counter-clockwise direction, 

as many times as there are unstable open-loop poles of KG'. For most practi- 
cal problems involving flexible structures, the only poles of G on the imagi- 
nary axis will be at the origin, and there will be an excess of poles over 

zeros, implying that limK(s)G(s)=0. In this case the Nyquist D-contour reduces 

s'*- 


1 Two other unstructured uncertainty representations are the additive and sub- 
tracted representations, which are discussed in Ref. [LEI]. These, however, 
are equivalent to the multiplicative and divided disturbances respectively 
and will not be discussed here. 

2 Another method for bounding the uncertainties is the cone-bounded perturba- 
tion [D04] . This essentially corresponds to a multidimensional scaling of 
the problem and is again equivalent to Eq. (4.2.8). 
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to the imaginary axis, with the possibility of an identation about the origin, 
and sefl R can be replaced by jw. Under the assumptions that G and G' have the 

same number of open-loop unstable poles, and that KG' and KG have identical 
open-loop poles on the jw-axis we have the following two results: 

Theorem 4.1: (Multiplicative Disturbance) The closed loop system will remain 
stable for all allowable perturbations iff: 

“o[KG(s)(I+KG(s))‘ 1 ] < l/l m (s) -¥sEn R (4.2.9a) 

or equivalently if KG is invertible: 

jL[I+(KG(s))” 1 ] > l ffl (s)*¥sea R (4.2.9b) 

Theorem 4.2: (divided disturbance) The closed loop system will remain stable 
for all allowable perturbations iff both a) and b) are true: 

a) L(s) has no zero or strictly negative real eigenvalues for any seft^ 

b) ZT[ (I+KG) -1 ] < l/l d (s) VseJ^ (4.2.10a) 

or equivalently: 

il[I+KG] > l d (s) -¥sen R (4.2.10b) 

A proof of the above results is given in Appendix B. 

It should be noted that these are not conservative results, given the error 
characterization of Eq. (4.2.8). In fact, if the above conditions are not 
met, there exists a perturbation, A ffi (s) or A d (s), whose maximum singular value 
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lies below l ffl (s) or l^Cs) respectively, which will destabilize the system. 

The unstructured uncertainty can however be unduly conservative in its charac- 
terization of error. If the error is known to have some structure it is quite 
possible that the particular A m (s) or A^(s) required to cause instability can- 
not occur. 

One major difference between the above two results is the additional 
requirement, in the case of a divided uncertainty, that L(s) have no zero or 
strictly negative real eigenvalues for ssfi R . This is needed to insure that 

the function [I+eA^] * remains continuous as t varies from zero to one, and 

can be interpreted as indicating that Theorem 4.2 cannot guarantee robustness 
in the case where phase is completely arbitrary. See Appendix B and 
Ref. [LEI] for further details. This places a limit on the situations in 
which a divided disturbance can be used to study robustness. For most practi- 
cal problems, however, the phase will be known to within ±180°, and both Theo- 
rems 4.1 and 4,2 are applicable. In this case the system can be made robust 
to multiplicative uncertainties by maintaining a large inverse return differ- 
ence, or equivalently a small loop gain, while the same system can be made 
robust to divided uncertainties by maintaining a large return difference, or 
equivalently a large loop gain. The two representations therefore imply oppo- 
site requirements for achieving robustness, and together imply that a system 
can be made robust either by maintaining very high loop gains, or by maintain- 
ing very low loop gains, while it will be most sensitive to errors in the 
region of gain cross-over. This is a familiar result from classical control, 
corresponding to the fact that perturbations of the Nyquist plot far away from 
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the -1 point will not affect stability, and the unstructured uncertainties 
provide its MIMO generalization. One final note is that Eq. (4.2.5b) indi- 
cates that the return difference and inverse return difference are not inde- 
pendent, implying that a system cannot be made robust to both multiplicative 
and divided disturbances acting simultaneously at any fixed frequency. 3 

In applying the unstructured uncertainties to flexible structures, which 
will have an infinite number of possibly lightly damped and poorly modeled 
modes, the divided disturbance is not particularily useful. The reason for 
this is that if some set of modes at "low frequencies" are robustly controlled 
by maintaining a high loop gain in that frequency range, there there will 
necessarily exist a further set of inodes at some "intermediate frequency" 
which will lie near cross-over. If the low frequency modes are inaccurately 
modeled, then the intermediate frequency modes will also be inaccurately 
modeled, and the system will not be robust to errors in these modes. The 
unstructured uncertainty, when applied to lightly damped flexible structures 
therefore implies that the loop gain must be kept below some level, whenever 
poorly modeled frequencies are present. To gain an appreciation for this 
limit, consider a very simple SISO flexible model. 


3 The divided disturbance is often considered most appropriate in dealing with 
low frequency errors, while the multiplicative disturbance is most appropri- 
ate when dealing with high frequency errors. This is because the divided 
disturbance indicates high loop gain, typically true at low frequencies, 
while the multiplicative disturbance indicates low loop gain, typically true 
at high frequencies. It should be emphasized that both are unstructured, 
with the only additional structure available in the divided disturbance 
being due to the fact that phase cannot be completely arbitrary. 
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SiaiEk °f i M ultiplicative aatructare.d Uncertainty Approach 

Sine, the unstructured uncertainty approaches are generalisations of c las- 
stcal SISO theory, they nil! work for more CMpllc « ed systems only 

worh fox a simple SISO system. With this in mind consider a system consisting 
of a single, underdamped, flexible mode whose gain and damping ratio are 

taown. hut whose frequency is Known only within given bounds. Such a system 
might be described as follows: 


2 


G’(s) =- T 


s*+2?f ' s +f ' 2 


f 2 

G(s) = 2 

s 2 +2£f s+f 2 


(4.2.11) 


let Af = (f . f o)/fo 


then A (s) = 

® s 2 + 2Cf' s + f' 2 

M m (jw)i = lAfiw 

(f -w 2 ) 2 + (2£f')V 


(4.2.12) 


(4.2.13) 


This ts the l m ( w ) for 
actually only known to be 


one particular error in frequency (Af), but Af is 
within certain bounds. The correct ljw) would then 


take the worst case of the above bound for 
find l a ( u) f or a set Qf bQunds Qn Af 

Eq. (4.2.12) for a sufficient number of Af' 
the worst case. 


every possible Af. An algorithm to 
for every «, calculate 

s between the bounds, and then plot 


Fig. 4.2a illustrates the result of such an algorithm. 
* 1 rad/sec, { = .01 and Af is allowed to vary between 


In this case 
-.10 and +.10. 
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201og 10 (l/l m (w)) is plotted, taking the worst case of 200 different Af's bet- 
ween -.10 and +.10 at 400 frequency points between .75 rad/sec and 1.25 rad/ 
sec. 

A computationally simpler approximation would be to consider only the two 
worst cases of Af, (Af * -.10 and Af * +.10). The result of an algorithm that 
did this for the same number of frequency points is plotted in Fig. 4.2b. The 
plots are similar, though the shapes vary slightly between the limits on Af. 
The second approach, might provide guidelines for a first cut design, though 
it would not strictly guarantee stability by Theorem 4.1. 

The function (1/1 («)) places strict limitations on the bandwidth of the 
m 

system, since the closed-loop gain must fall below it. This implies that the 
bandwidth is limited by the first time that (1/1 («)) falls below 0 db. This 

ID 

occurs when l m (w) first rises above 1, or when the errors in frequency res- 
ponse first rise above the nominal frequency response. It would be expected 
that this would occur at some frequency for any realistic plant but will it 
occur near a given uncertain natural frequency? In the cases plotted in Figs . 
4.2a and 4.2b (l/l m (w)) actually dropped to about -20 db. Theorem 4.1 would 

therefore imply that a SISO feedback system with with a 1% damping ratio and 
10% frequency error would be required to have a loop gain that fell below 
-20 db at that frequency, sharply limiting bandwidth and therefore perfor- 
mance. 

To get a rough idea of the extent of the limits imposed on such a system by 
Theorem 4.1, consider an approximation of Eq. (4.2.13). Assume that £ and Af 
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are both small as compared to 1, and that the function reaches a maximum at 
«_= f ' . Then: 

|A ( j«)l max a Af/C (4.2.14) 

This indicates that relative frequency errors on the order of the damping 
ratio will limit the system closed-loop bandwidth. This is a relatively 
strict limitation, since if the damping ratio in a particular mode is 1%, a 
reasonable value for flexible space structures, a 1% frequency error might 
destabilize a system whose bandwidth included that mode. 

This limit seems overly restrictive, but to gain a more concrete apprecia- 
tion of the conservativeness of the unstructured uncertainty for this particu- 
lar system consider a simple constant gain controller in the feedback loop. 

The closed-loop characteristic equation of the system is: 

s 2 + 2Cf's + (l+k)f' 2 = 0 (4.2.15) 

This is clearly stable for all k > -1, for all C > 0 and for all f'. The 
bandwidth of the closed- loop system can therefore be increased without bound, 
for arbitrary error in frequency and arbitrary non-zero damping ratios, in 
conflict with the requirements implied by Theorem 4.1. The reason that the 
unstructured uncertainty is so conservative, even for this very simple, SISO 
system, is that it doesn't take into account any information on the phase of 
the system. Essentially it defines limits on the gain error, while allowing 
completely arbitrary phase error. For this system though, the phase error is 
uniquely determined by the gain error, so that the exact phase shifts that 
would be necessary to destabilize the system for a given gain shift can never 
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occur. The unstructured uncertainty is therefore not a good characterization 
of error due to uncertain frequencies in the modeled modes of a lightly 
damped, flexible structure. Next apply Theorem 4.2 to the same problem. 


Application of Theorem 4.2 to the Example Problem 

In applying Theorem 4.2, the first step is to ascertain that L(s) is never 
zero, or strictly negative on seft^, or equivalently that the phase of L(ju) 

never reaches ±180°. For the example: 

f' 2 [(f 2 -w 2 )(f’ 2 -w 2 )+4c 2 f fV+2£w[f n (f' 2 -w 2 )-f7f 2 -« 2 )]j] 

L(jw) = - (4.2.16) 

w 2 [ (f ' 2 -w 2 ) 2 +4£ 2 f ' 2 w 2 ) 

L(jw) is equal to zero or strictly negative only when the imaginary part of 
eq. (4.2.16) is identically equal to zero. This in turn is true only when £=0 
or when w^-f^f 1 , so unless C=0, there is no real frequency at which L(jw) is 

either zero or negative real, indicating that Theorem 4.2 is applicable . k For 
the example the divided disturbance is given as follows: 


M*) = 


! +2{f 0 s+f 0 


where Af=(f Q -f ' )/f ' 


(4.2.17) 


|Aj(jw)| will reach a maximum near w=fQ, so: | A^ ( jw) l m ay = Af/C 


Theorem 4.2 therefore implies that whenever the frequency error is on the 
order of the damping ratio, loop gain at that frequency must lie above Odb, or 


* Actually, for small £ and large Af the phase of L(jw) asymptotically reaches 
±180* , but it is never exactly equal to ±180*. 
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equivalently, that sufficiently high gain will stabilize the system. Theorems 
4.1 and 4.2 together indicate that a feedback control for this example can be 
made robust by either introducing very high gains, or very low gains, while in 
the case where Af>C there will exist a zone of intermediate gains which will 
not be robust, as was pointed out in the previous Section. This is a conser- 
vative result since it has been ascertained that the above example is robust 
to arbitrary frequency errors for any positive feedback gain. 

In conclusion, errors in the frequency of a lightly damped oscillator pro- 
vide a highly structured parameter variation and the unstructured uncertain- 
ties of Doyle and Lehtomaki provide an overly conservative characterization of 
these errors. This further implies that due to their conservativeness 
unstructured uncertainties are not particularly useful for determining the 
robustness of flexible systems with poorly modeled frequencies. An implica- 
tion of this observation will be dicussed further in the context of LTR design 
methods (Section 4.2.4). 


4.2.2 ASYMPTOTIC LOOP TRANSFER RECOVERY (LTR) CONTROL DESIGN 

LTR is an appealing approach to control design for two basic reasons. The 
first is that it can be used to recover the gain and phase margins of a Linear 
Quadratric Regulator (LQR) at the plant input, or of a Kalman Bucy Filter 
(KBF) at the plant output. This is sometimes called robustness recovery [D02] 
or sensitivity minimization [KW1], but these are somewhat misleading terms 
since they suggest that the loop recovered system will have the same robust- 
ness properties as the full state feedback LQ regulator or the KB filter. 

This is not necessarily true as will be indicated more specifically later. 
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The second appeal of the LTR approach is its applicability to loop shaping as 
a more general control design method. The idea of loop shaping is most 
clearly spelled out in Ref. [D01], where it is pointed out that a number of 
system properties, including performance, sensitivity, noise rejection, 
"robustness" and control effort, depend on the system loop shape (maximum and 
minimum singular values of the system transfer function matrix). The LTR loop 
shaping approach involves first designing an LQR loop or a KBF loop with 
desired characteristics, and then recovering that loop shape by LTR methods. 
Again there is an assumption in the loop shaping approach that all systems 
with the same loop shape will behave identically. This is true for the nomi- 
nal plant, but when the plant parameters are perturbed, the two loops may no 
longer be identical and the system response (and stability) may vary consider- 
ably. In a later Section, a number of examples will be given of systems with 
identical loop shape, but different robustness characteristics. 

For simplicity, the input loop recovery procedure of Doyle and Stein 
[D01,2] will be considered, noting that output loop recovery [KW1] is simply 
its dual. After describing the procedure it will be interpreted for SISO sys- 
tems, and for 1-input, m-output systems. This interpretation will lead to a 
polynomial (pole/zero cancellation) approach to loop recovery. Methods for 
using the extra degrees of freedom available in the loop recovery procedure to 
find an "optimal" loop recovered system will be considered in Section 4.4. 


Recovery 


The LQR loop has excellent robustness properties. In particular an LQ 
regulator with diagonal input weighting matrix R is guaranteed to have at 
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least an infinite gain margin, a gain reduction margin of 1/2 and a phase mar- 
gin of ±60°, simultaneously in all the feedback loops [LEI, 2]. The idea of 
input loop recovery is to design an LQ regulator with a desired loop shape and 
then use an asymptotic procedure to design a KB filter that recovers that loop 
shape. Since gain and phase margins are a function of loop shape, the output 
feedback system, using a KB filter to estimate the plant states, will have the 
identical gain and phase margins as the LQR loop. For this reason, input loop 
recovery is sometimes called robustness recovery. This is the point of view 
taken in Ref. [D02] . Another advantage of the loop recovery approach is that 
it is relatively simple to specify the LQR loop shape at low frequencies by 
carefully choosing the weightings in the LQR cost functional [HA1] . Loop 
Transfer Recovery (LTR) can then be used to achieve the identical loop shape 
for an output feedback system. This is the loop shaping control approach, and 
is the point of view taken in Ref. [D01]. 

To see how loop transfer recovery works, consider a full-state feedback 
control law, u=-Kx. The input loop transfer function for this system will be 

K(sl-A) ^B. Now consider a state estimator of the form, 

$ = A$ + Bu +G(y - Cx) (4.2.18) 

and a control law based on the state estimate 5, u=-K$. In this case the com- 
pensator transfer function is, 

K(s) = K(sl - A + BK + GC) -1 G (4.2.19) 

and the input loop transfer function for the output feedback system is: 

K(s)G(s) * K(sl - A + BK +GC) _1 GC(sI - A) -1 B (4.2.20) 

To get loop recovery choose an estimator gain matrix G, such that, 

K(sI-A+BK+GC)" 1 GC(sI-A)' 1 B K(sI-A) -1 B (4.2.21) 
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In Refs. [D01,2] it is indicated that any G matrix that asymptotically 
satisfies 

G -► qBW (4.2.22) 

for some arbitrary symmetric, positive definite matrix W, will recover the 
loop shape, q is a large scalar parameter which is increased to achieve 
asymptotic loop recovery. To illustrate this result, note that the transfer 
function between y and y in the state estimator can be written as follows, 

$ = (I + ClG^ClGy (4.2.23) 

where i denotes (si - A + BK) This leads to the following alternate repre- 
sentation of K(s), 

K(s) = K[7 - ?(I + ciG)' 1 cl]G (4.2.24) 

which after some further manipulation can be rewritten as: 

K(s) = k7g(I + C?G) -1 (4.2.25) 

Now let G *+ qBW which means that, 

K(s) - k7b(C?B)' 1 . (4.2.26) 

Using the identity, $B = §B(I + K#B) *, where i = (si - A) 1 Eq. (4.2.26) 
becomes , 

K(s) KfB(CfB)' 1 (4.2.27) 

and, 

K(s)G(s) - K*B(C#B) -1 C*B = K$B (4.2.28) 

This shows explicitly how loop recovery inverts the plant dynamics from the 
left and replaces them by the full-state feedback dynamics. 

So far there has been no requirement that K be derived from an LQR problem, 
nor that the estimator be a KB-filter. In fact any appropriate method (e.g. , 


48 


pole placement) could be used to select the control gains K, and any estimator 
(e.g., an observer) which both asymptotically satisfies Eq. (4.2.22) and sta- 
bilizes the plant, could be used to calculate the asymptotic estimator gains 
G. However, the choice of an LQR approach for finding K has the advantage 
that the' resulting loop will have desirable properties, and the use of a KB- 
filter to recover that loop shape has the advantage of guaranteeing a stable 
plant at every stage in the asymptotic loop recovery procedure. 

To use a KB filter for the LTR, the first step involves appending addi- 
tional columns to the B matrix until the plant is square. The only require- 
ment on these columns is that the resulting plant be minimum phase, and within 
this constraint the extra columns can be viewed as free design variables. 

Once the system is "squared up," a KB filter is found, where the measurement 
noise covariance N is an arbitrary positive definite matrix, and the state 
noise covariance has the following special form, 

M = M q + q 2 BWB T (4.2.29) 

where is some nominal noise covariance, q is the parameter which will be 

increased to achieve recovery and W is an arbitrary positive definite matrix. 
Note that the B matrix in Eq. (4.2.29) includes the appended columns, so the 
T 

rank of BWB is equal to the number of plant outputs. 

To show that this KB filter will asymptotically satisfy Eq. (4.2.22), con- 
sider the Riccati equation: 

AP + PA T - PC T N -1 CP + M q + q 2 BWB T = 0 (4.2.30) 

2 

and divide by q : 
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(4.2.31) 


A(P/q 2 ) + (P/q 2 )A T - q 2 (G/q 2 )N(G/q 2 ) T + (M Q /q 2 ) + BWB T = 0 
Now let q «•, under the assumption that the plant be minimum phase, to get: 

-(G/q)N(G/q) T + BWB T -*■ 0 (4.2.32) 

The KBF estimator gains will therefore have the following asymptotic property: 

G qBW*TN~* (4.2.32) 

T T 

Where T is an orthonormal matrix (ie. TT = T T s I). This of course satis- 
fies Eq. (4.2.22) and the KB filter will asymptotically recover the loop shape 
as q is made to approach infinity. 

One quick note before continuing to an interpretation of loop recovery for 
SISO systems, is that the estimator gain matrix G may become very large before 
loop recovery is achieved. This may cause some computational problems and may 
not provide the best solution for a robust compensator, as will be discussed 
further in Section 4.2.3. A related problem concerns the convergence of esti- 
mator functional gains. As G increases in size, the performance of the esti- 
mator increases, implying that the number of modes necessary to achieve con- 
vergence of the estimator functional gains also increases. Asymptotic LTR 
therefore implies a large system model if the estimator is expected to con- 
verge to the optimal infinite dimensional estimator, as measured by functional 
gains. This issue still needs to be looked into. Next the SISO case will be 
examimed to gain a more intuitive understanding of the LTR procedure. 


SISO Interpretation of Loop Recovery 

The SISO loop recovery problem is as follows. Given a SISO 
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G(s) * n(s)/d(s), find a proper (or strictly proper) K(s) = N(s)/D(s) such 
that K(s)G(s) = G(s)K(s) * 3f(s)/d(s) for some specified Y(s). The obvious 
solution is to let N(s) * 2f(s) and D(s) = n(s)D^(s), where D^(s) consists of 

enough poles to make K(s) proper (or strictly proper), and these poles are 
placed far enough into the l.h.p. that they do not significantly affect the 
loop shape in the design region. This interpretation also clarifies the mini- 
mum phase condition, since if G(s) is non-minimum phase, this procedure will 

result in a pole/zero cancellation in the r.h.p., that will destabilize the 

closed loop system. 

The asymptotic LTR procedure does not calculate an exact pole/zero cancel- 
lation, but achieves this cancellation asymptotically. It is clear, however, 
that it is possible to achieve loop recovery by specifying a compensator that 
does achieve exact cancellation, without going through an asymptotic process. 
In both cases the cancellation will be close to exact before loop recovery is 
achieved, so it will be assumed that it is exact. 

This interpretation indicates the essential difference between the true LQR 
loop and the LTR loop. While both loops may look identical over any given 

frequency range, the LTR loop will contain a number of hidden pole/zero can- 

cellations. These cancellations do not show up for the nominal plant, but as 
soon as the plant changes, they will no longer be exact and the LQR and LTR 
loop shapes may be considerably different. This is especially evident for 
lightly damped systems, since the plant zeros in this case will lie close to 
the imaginary axis. Very small errors in the plant zeros will therefore pro- 
duce very large errors in the LTR loop shape. It is because of this property 
that an LTR design will in general not have the same robustness character is - 
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tics as the corresponding LQR design. 


Some Comments on the MIMO Case 

LTR is a more interesting and less obvious design procedure in the MIMO 
case, but the same basic properties carry through. In this case the designer 
would choose the following compensator, 

K(s) = K§B(C§B)" 1 /D 1 (s) (4.2.34) 

where D^(s) would have the same interpretation as the SISO case. Eq. (4.2.34) 

can be rewritten by dividing out the common factor det[sl - A] to get, 

K(s) = (Kadj(sI-A)B] [Cadj(sI-A)B]‘ 1 /D 1 (s) (4.2.35) 

where Kadj(sI-A)B is the MIMO generalization of Jf(s) and Cadj(sI-A)B is the 
MIMO generalization of n(s). 

For the special case of 1-input and m-outputs [Cadj (sI-A)B] will be a 
square matrix where only the first column is fixed and the next m-1 columns 
are at the designers discretion (within the constraint of a non-minimum phase 
plant) . Thinking again in terms of the SISO problem, this is somewhat like 
having an n(s) at the designer's discretion, corresponding to some freedom in 
the placement of compensator poles. Once the poles are picked, however, loop 
recovery will specify the compensator zeros. Keeping this in mind an alge- 
braic, direct pole/zero cancellation, design of a LTR compensator for a 
1-input, m-output plant will be considered. 
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4.2.3 ALGEBRAIC LOOP TRANSFER RECOVERY 


As indicated in the previous Section, loop recovery can be achieved for a 
SISO system by defining a compensator that exactly cancels the plant zeros and 
replaces them by a set of "optimal" zeros corresponding to the LQR loop. For 
the case of a 1-input, m-output plant, a similar procedure can be carried out, 
but now the designer has a set of free design variables. This freedom can be 
used in a variety of ways. Some possibilities include choosing a minimum 
order loop recovery compensator or specifying the location of the compensator 
poles. For the purposes of this discussion a minimal order compensator for 
which the compensator poles can be chosen arbitrarily will be considered, and 
an algebraic procedure to find the compensator numerator polynomials, or equi- 
valently the compensator zeros will be used to achieve a given input loop 
shape. It should be emphasized that this is not a requirement of the alge- 
braic LTR approach. Lower order compensators might be found by removing the 
freedom to specify pole location, while higher order compensators would allow 
the designer a number of free parameters which could be adjusted so as to meet 
some other criteria, possibly via an optimization process. The choice consid- 
ered for this discussion, however, allows one to examine the important effect 
of compensator pole location, without adding excess variables to the formula- 
tion. 

A 1-input, m-output plant and compensator will have the following special 
form: 
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G(s) = 


gjCs)" 


’n 1 (s)" 

g 2 (s) 


n 2 (s) 

• 

S 

• 



n (s) 

l no J 


/d(s) 


K(s) = [k^s) k 2 (s) ... k m (s)] = [NjCs) N 2 (s) ... N m (s)]/D(s) 


The input loop transfer function will be: 


(4.2.36a) 


(4.2.36b) 


m 

K(s)G(s) = 2n.(s)N.(s)/d(s)D(s) (4.2.37) 

i=l 

To achieve loop recovery it is necessary that: 
m 

K(s)G(s) = £ n (s)N (s)/d(s)D(s) = T(s)/d(s)D. (s) = I(s)/d(s) (4.2.38) 

i=l 1 1 

where 2f(s) is the numerator polynomial for the LQR input loop. D^(s) is again 

a set of poles required to maintain a proper (or strictly proper) compensator, 
where these poles are kept far enough out in the l.h.p. that they don’t signi- 
ficantly affect the loop shape in the design region. Then D(s)=D c (s)D^(s) , 

where D c (s) must be cancelled out by the input loop numerator. Once the com- 
pensator poles are selected, the following equation must be solved to find the 
compensator numerator polynomials: 

m 

2 n i (s)N.(s) = »(s)D c (s) (4.2.39) 

i=l 

Next consider the minimal order compensator that can satisfy Eq. (4.2.39). 
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Let xiy be deg(T(s)), be the degCD^Cs)), be deg(D^(s)) and n N be 
deg(N^(s)). The degree of the left side of Eq. (4.2.39) will be + n fl , 
while the degree of the right side of Eq. (4.2.39) will be n Y + n , so the 
first requirement on n^ is that: 

“N + n n = n Y + n c (4.2.40) 

The second requirement on n^ is that there must exist enough degrees of free- 
dom to satisfy Eq. (4.2.39) for arbitrary y(s)D (s). The number of free par- 

c 

ameters is m(n^+l), while the number of polynomial coefficients that must be 
satisfied is n„+n +1, this implies that: 

4 C 

m(n N +l) = n y + n c + 1 (4.2.41) 

Combining Eqs. (4.2.41) and (4.2.42) gives the following result for n^: 

= (n n /(m-l)) - 1 (4.2.42) 

For a strictly proper compensator, the order of the compensator is 
n + n 1 = n /(m-1). So a 1-input, 2-output plant has a strictly proper loop 

C i n 

recovery compensator of the same order as the degree of the plant numerator 
polynomials, while a 1-input, 3-output plant has a loop recovery compensator 
of half that order, etc.. This approach therefore leads to lower order com- 
pensators than the asymptotic, observer based approach described in the previ- 
ous Section. In fact, as m increases, the order of the compensator can become 
very small, while still achieving loop recovery. 

Once the order of the compensator and the location of the compensator poles 
are chosen, the compensator numerator polynomials are found by the solution of 
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a set of m(n n +l) linear equations. Let the ith plant numerator polynomial be 

written as follows: 
n 


n.(s) = n. s + ... + n. ,s + n. . 

1 ' i,n i,l i,0 

n 

and similarily let the ith compensator numerator polynomial be: 


n 


N. 


N.(s) = N. s + ... +N. , s + N . A 
l i.njq i,l i,0 


Then: 


(4.2.43a) 


(4.2.43b) 
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(4.2.44\ 


n +n 
n nn 


or symbolically, Nn = d, where d. is the ith coefficient of 2T(s)D (s). The 

i c 

above algorithm is easily implemented on digital computer whenever a unique 
solution to Eq. (4.2.39) exists, or equivalently whenever the matrix N is 
non- singular. 
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4.2.4 ROBUSTNESS OF LTR DESIGNS 


The purpose of this Section is to present a discussion of particular 
aspects of LTR designs which relate to robustness. 

One observation concerns the relationship between loop shape and robustness 
implied by Theorems 4.1 and 4.2. Specifically, if model errors are accurately 
described by unstructured uncertainties, then LTR methods can be used to 
ensure that the conditions of Theorems 4.1 and 4.2 are met, thereby guarantee- 
ing robust stability. This constitutes the ideal case for applying Loop 
Transfer Recovery, since specifying constraints on the loop shape is a neces- 
sary and sufficient condition for robustness. On the other hand, some model 
errors are not well described by unstructured uncertainties. In particular 
frequency errors in the model of a lightly damped flexible structure fall into 
this category (see Section 4.2.1), as do most cases of strict parameter varia- 
tions. In this case Theorems 4.1 and 4.2 give sufficient but not necessary 
conditions for robust stability, and a design based upon constraining loop 
shape on the basis of unstructured uncertainties will be conservative, result- 
ing in low performance. In the case of lightly damped structures the implied 
constraints on loop shape would require that uncertain flexible modes not be 
actively controlled (again see Section 4.2.1). Therefore, the robustness of 
designs which actively control flexible modes cannot be analyzed by examining 
loop shape alone. In fact, for a given loop shape, a whole family of LTR 
designs may exist, where the robustness of these designs varies. In conclu- 
sion, the application of LTR methods to lightly damped structures with uncer- 
tain frequencies cannot be blindly approached in terms of conditions implied 
on the loop shape. A major emphasis of the current study has been to identify 


57 



conditions which are important. These results will be illustrated by the 
examples presented in Section 4.2.5. 

As noted in Section 4.2.2, one variable among a family of LTR designs with 
the same loop shape is the location of pole/zero cancellations. Whenever 
these cancellations occur near the imaginary axis a small perturbation in the 
plant, resulting in a small shift of the plant zeros, will cause a very large 
variation in the open- loop frequency response. It would therefore be reason- 
able to expect that an LTR design which tries to avoid pole/zero cancellations 
near the imaginary axis will be more robust than one that does not. One way 
to reduce the number of plant zeros, and therefore the number of cancella- 
tions, is to add measurements. A second order flexible system, with m-inputs 
and m-displacement outputs, can have at most (n-2m) transmission zeros [EDI] 1 . 
Remembering that in the asymptotic LTR approach, each additional measurement 
entails an additional artificial input, this implies that every measurement 
removes a pair of transmission zeros. In the version of the algebraic 
approach used at this point, additional measurements are used to reduce the 
compensator order, so it is unclear whether any robustness improvement might 
result. 

A method for shifting the cancellations away from the imaginary axis in the 
asymptotic approach would be to use the extra degrees of freedom available in 
the choice of artificial inputs to shift the plant transmission zeros to the 
left. The best way to do this is still under investigation. This idea is 
also the motivation behind the algebraic LTR approach, which allows the desig- 

1 For velocity outputs the maximum number of finite transmission zeros is 
(n-2m+l) 
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ner to place the pole/zero cancellations as far into the left half plane as 
desired. The general question of how to choose the most robust LTR compensa- 
tor for a given LQR loop shape is certainly not closed, and one method which 
we will continue to study is the optimal sensitivity approach described in 
Sec. 4.4. 

One final comment on the robustness of LTR designs is a short discussion of 
Ref. [SHI], which deals specifically with the robustness problems associated 
with any compensator design that incorporates very high estimator gains. 
Asymptotic LTR designs can result in extremely high (on the order of 10* or 
higher) gains in the estimator Kalman gain matrix G. Two difficulties may 
arise. The first involves the case where a plant perturbation changes the 
structure of the plant transmission zeros at infinity. This is a perturbation 
that affects the excess of poles over finite transmission zeros 2 . An example 
in which this occurs for a realistic flexible system is as follows: Consider 

the system illustrated in Fig. 4.3, consisting of two disks (of inertia 1), 
stacked one on top of the other, connected by a rod (of length 1) which acts 
as a rotational spring (of spring constant 1/2). The input is a torque 
applied to the lower disk. The nominal output is the rotation of the upper 
disk. The nominal transfer function is: 

1/2 

G(s) = s 2 (s*+l) (4.2.45) 


2 Note that this is not the same as the finite transmission zeros added when 
neglected dynamics are appended to the design plant. In this case the num- 
ber of poles increases along with the number of transmission zeros, indicat- 
ing that the excess does not change, and the structure of the zeros at 
infinity therefore does not change. 
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An LQR design which minimizes the mechanical energy, plus the rotation of the 
lower disk squared, plus the control input squared, results in the following 
control gain vector (in modal coordinates): 


K = [ 1.0000 2.0555 1.6124 0.6325 ] 


(4.2.46) 


An asymptotic LTR design can be achieved by setting the process noise covari- 

12 T 

ance equal to (1x10 )BB , resulting in the following Kalman Filter gains: 


G = 


98,255.0’ 

1 , 000 , 000.0 

98,116.0 

995,170.0 


(4.2.47) 


This design is relatively robust with respect to frequency variations (with- 
stands ±30% errors in the flexible frequency), but consider an error corres- 
ponding to a manufacturing defect which places the rotation sensor at a point 
t below the top disk. The true transfer function would then become: 


r . M _ es*+l/2 
G (S) " s 2 (s 2 +l) 


(4.2.48) 


Therefore, while the nominal plant has no finite zeros, the true plant has a 
pair at ±/-l/2e, satisfying the conditions referred to in Ref. [SHI]. The 
asymptotic LTR design defined by Eqs. (4.2.46) and (4.2.47) will be unstable 
for e>0.04, 3 which produces finite zeros at ±3.54j or lower on the imaginary 
axis. This might be a troublesome robustness problem in some applications, 
but it can occur only under very special circumstances. For any square flexi- 


3 It will also be unstable for e<-0.01, but this corresponds to a measurement 
location above the disk and is not physically feasible. 
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ble structure with displacement outputs, the maximum number of finite tran- 
smission zeros is (n-2m) [EDI]. If the nominal plant has this maximum number 
of transmission zeros, then no perturbation in plant parameters can result in 
the conditions required by Ref. [SHI] for robustness problems. 

In conclusion the robustness problems associated with LTR control exposed 
in Ref. [SHI] depend on a nominal plant which is chosen such that there exists 
a perturbation in plant parameters which will increase the number of finite 
transmission zeros by at least 2. This possibility can always be ruled out by 
ensuring that the nominal (design) plant has the maximum possible number of 
finite transmission zeros, which will be the case in most practical problems. 
In particular the antenna model studied in this report always has the maximum 
number of transmission zeros and is not affected by the robustness problems 
exposed in Ref. [SHI]. The extension of these ideas to infinite-dimensional 
systems is not considered in Ref. [SHI] and has not been addressed in the pre- 
sent work. This would, however, provide an interesting topic for future exa- 
mination. 

The second observation in [SHI] is that, due to the very high gains 
involved, a small error in the realization of the compensator could destabi- 
lize the system. The authors suggest that a compensator be realized directly 

from the transfer function K(sI-A+GC+BK) *G. This issue also has not been 
considered in detail, but it will be examined later. 
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4.2.5 EXAMPLES OF THE LTR APPROACH 


Antenna Model 

The methods developed in the previous Section are applied to the wrap-rib 
antenna described in Chapter 2. The system is modeled by a component mode 
method, ie. mode shapes of the beams and the mesh Sections are found sepa- 
rately, and the "modal" representations are then combined to form a model for 
the 90° sector. A subset of six modes are selected from this very large model 
to constitute a design model for one quadrant of the antenna Frequencies, 
modal input influence coefficients, and modal output influence coefficients 
for the six modes included in the design model are listed in Table 4.2.1. 

Damping is assumed to be visco-elastic (ie. the damping matrix is propor- 
tional to the stiffness matrix) 1 , but the model allows different damping coef- 
ficients 2 for the beam and the mesh portions of the sector model. This 
implies that damping is added at the component mode stage, and furthermore 
that the damping matrix, after combining the beam and mesh component mode 
models and calculating the overall modal representation, will include off- 
diagonal terms. The final antenna model, in a second-order modal representa- 
tion, will have the following form: 

V + D0 + A 2 # = Bu, y ■ C$ (4.2.49) 

where A 2 is the diagonal matrix of the natural frequencies squared, while D is 


1 Visco-elastic damping implies that the damping ratio of a mode increases 
proportionally with the frequency of that mode. Higher frequency modes will 
therefore be more highly damped. This is a realistic situation for most 
materials, particularily in the higher modes. 

2 The damping coefficient is a scalar number which multiplies the stiffness 
matrix, resulting in the damping matrix. 
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symmetric, positive-semidefinite, but is diagonal only in the case when the 
damping coefficients in both the beam and the mesh are identical. The parti- 
cular model studied here has higher damping in the mesh, and the damping 
matrix is therefore not diagonal, though the off-diagonal terms are in general 
an order of magnitude smaller than the diagonal terms. The entire damping 
matrix is listed in Table 4.2.1. 

Table 4.2.1 Antenna Modal Data 


Mode # 

1 

2 

3 

4 

5 

6 

w 

0 . 

6.95 

18.95 

36.57 

51.89 

55.32 

C 

0 . 

.0105 

.0302 

.0589 

.0519 

.0537 

B 

o 

.093 

-.122 

-.285 

-.511 

-.120 

-.111 

C 1 

.093 

-.122 

-.285 

- .5H 

-.120 

-.111 

C 2 

.0 

.243 

.183 

.617 

.201 

.115 

C 3 

.0 

.002 

.183 

.341 

.099 

.052 

C 4 

.0 

.001 

.129 

.241 

.070 

.037 

C 5 

.0 

.000 

-.019 

.068 

-.027 

.020 

C 6 

.0 

.172 

.129 

.437 

.142 

.081 


Damping Matrix 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.146 

-0.030 

-0.050 

-0.016 

-0.009 

0.0 

-0.030 

1.144 

-0.151 

-0.096 

-0.049 

0.0 

-0.050 

-0.151 

4.308 

-0.009 

-0.095 

0.0 

-0.016 

-0.096 

-0.009 

5.385 

0.099 

0.0 

-0.009 

-0.049 

-0.095 

0.099 

5.940 


Simpler models, including series of connected masses and springs, as well 
as a beam/hub model have also been studied, but the difficulty in achieving 
robust control designs has been much greater for the antenna sector model. It 
is hypothesized that this is due largely to the addition of mesh modes to the 
beam/hub model, resulting in a densely packed set of frequencies, and some 
very slightly controllable/observable modes. From the point of view of con- 
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trol design, the important aspect of this particular model is that it is 
difficult to control robustly, and therefore poses a challenging problem. 

The number of available measurements (or plant outputs) is not fixed at 
this time. There must be at least one measurement relative to inertial space, 
to insure that the rigid body mode is observable. This measurement is pro- 
vided by a rotation sensor colocated with the hub torque input. To effec- 
tively control the shape of the antenna, measurements on the antenna surface 
must also be taken. Possible locations for these measurements are labeled and 
numbered on Fig. 4.4, where all these measurements are relative to the central 
hub position. The modal output influence coefficients are listed in 
Table 4.2.1. In the compensator design stage, the measurements used will be 
referred to by the corresponding number found on Fig. 4.4. For the examples 
presented here, only one other measurement will be used. That measurement 
corresponds to position # 2, at the tip of the beam which lies perpendicular to 
the axis of rotation. 

One final note on the measurement locations is that all seven measurements 
are not necessarily independent for a given mathematical model. This means 
that, for a particular mathematical model studied, the value of one measure- 
ment may be directly related to that of another measurement. As an example, 
consider measurement locations 2 and 6. The complete mathematical model of 
the antenna sector includes symmetric beam modes (the two beams move together 
in the same direction) and asymmetric beam modes (the motion of one beam 
exactly opposes the motion of the other). With no mesh, the asymmetric modes 
are uncontrollable, and even with the mesh they are only very slightly con- 
trollable. These modes were therefore neglected in the mathematical design 
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model, since they do not have much affect on the control problem. This, 
however, implies that in terms of the mathematical design model, the two beams 
will always move together, and the measurements at locations 2 and 6 will be 
linearly dependent upon each other. This points to an important consideration 
in determining the independence of measurements. The fact that measurements 
are independent for a physical system or for some mathematical model does not 
imply that they are independent for every possible, reduced order, mathemati- 
cal model. Measurements 2 and 6 are clearly independent for the overall 
mathematical (and the physical) antenna model, but due to the exclusion of 
asymmetric modes, they are not independent for the reduced order mathematical 
design model studied here, due to the exclusion of asymmetric modes. In terms 
of this study, the implication is that adding measurement number 6 to number 2 
will not result in any further information on the state of the system. This 
is particularily important with respect to the algebraic design method which 
must assume independence of all the measurements used. 


LQR Design Problem 

The first step in a Loop Transfer Recovery (LTR) control design is the 
choice of a full -state feedback control law, or equivalently the matrix K of 
control gains. The matrix K can be found by a number of methods, including 
pole placement, but, perhaps the most advantageous is by solving a Linear 
Quadratic Regulator (LQR) problem. The statement of the LQR problem is as 
follows : 

Given x = Ax + Bu (4.2.50) 
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Minimize the following cost functional: 

(x^Qx + u T Ru)dt (4.2.51) 

0 

where Q is symmetric, positive semi-definite and R is symmetric, positive 
definite. The solution to the problem, for full state feedback is: 

u = -Kx (4.2.52) 

-1 T 

where K = R BP and P satisfies the following Riccati equation: 

PA + A T P - PBR -1 B T P + Q = 0 (4.2.53) 

The LQR approach has a number of positive features, including guaranteed 
gain and phase margins [LEI, 2] 3 , the ability to shape the loop gain at low 
frequencies [HA1] and the simplicity of the approach. The designer "simply" 
chooses the Q and R matrices and the computer does the rest . 

For the antenna problem, the input is a scalar, so without loss of general- 
ity R can be replaced by 1. The problem then reduces to finding a Q matrix in 
the following form: 

Q = q c Q 0 (4.2.54) 

where q c is a scalar parameter which is increased to increase system perfor- 
mance (and loop gain) , and Q q is a matrix which specifies the form of the cost 

weighting on the states. There are a number of bases for choosing Q. One 
might wish to achieve a particular loop shape, or a particular closed loop 

3 This constitutes a large part of the appeal of the LTR control design 
approach 
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pole configuration, or there might be some considerations in terms of a 
physically meaningful objective function to be minimized. For instance, one 
might want to minimize the mechanical energy in the system, or a pointing 
error, or perhaps the RMS error of a surface or signal. 

Three LQR designs will be presented for the antenna. The first minimizes 
the RMS surface error, since this provides a good meausure of antenna perfor- 
mance. The second minimizes the mechanical energy in the antenna, plus a term 
corresponding to the rigid body mode (this is necessary since there is no 
mechanical energy corresponding to rigid body rotation) . The third minimizes 
RMS error, plus the rate of change of RMS error. All three LQR designs remain 
stable for combined frequency errors of -99% or +100%. Any robustness diffi- 
culties in the LTR designs therefore are not due to similar difficulties in 
the corresponding LQR design. 


LQR Design it 1 


The first LQR design penalizes the RMS surface error of the antenna. This 
was chosen because the antenna performance is expected to be closely related 
to the surface error. The Q matrix which penalizes RMS error has the follow- 
ing form: 


Q 1 



0 

0 


(4.2.55) 


The elements of Q r are given in Table 4.2.2. Notice that the weightings on 

the rigid body mode and the first three flexible modes are approximately 
equal, while the weightings on the last two modes are higher by approximately 
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a factor of 5. Q f therefore approximately weights the position of the rigid 

body mode and the first three flexible modes equally, places a little more 
emphasis on the last two modes, and places no weighting on the modal veloci- 


ties. The resulting closed loop pole configuration, as q is increased from 

c 

.01 to 100, is illustrated in Fig. 4.5. An asterisk is placed by the closed 
loop poles for q = 100, and the loop gain for this case is illustrated in 

G 

Fig. 4.6. This loop gain indicates flexible modes both above and below Odb. 
This is a desirable, yet realistic situation, since if all modes were below 
Odb they would essentially be uncontrolled, while if an modes were above Odb, 
there would have to be some modes (of the true system) below Odb. 

Table 4.2.2 Matrix Q 

r 


“ 2166.16 

385.25 

-313.39 

-754.52 

-5980.06 

-4880 . 05* 

385.25 

1618.56 

-632.98 

42.80 

-5976.23 

-3724.72 

-313.39 

-632.98 

1598.49 

-771.90 

6803.32 

3031.96 

- 754.52 

42.80 

-771.90 

2539.28 

-5295.22 

1537.05 

-5980.06 

-5976.23 

6803.32 

-5295.22 

98137.61 

-2965.06 

-4880.05 

-3724.72 

3031.96 

1537.05 

-2965.06 

102461.54 


LQR Design #2 

The second LQR design penalizes the mechanical energy (kinetic + 
potential), plus a term corresponding to the rigid body mode. For a 

T 

mechanical system the kinetic energy is given by x Mx, where M is the mass 

T 

matrix, and the potential energy is given by x Kx, where K is the stiffness 
matrix. The appropriate Q matrix to weight the mechanical energy is 
therefore: 
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(4.2.56) 


K 0 

0 M 

If the rigid body coordinate is the first state variable, K will have a zero 
first row and first column, indicating that there is no potential energy cor- 
responding to rigid body motion. The resulting control law would therefore 
fail to control rigid body motion, so some weighting on the rigid body 
motion must be added to the K matrix. The new K matrix will be identical 
with the addition of a scalar parameter in the (1,1) position (q^ can be 

used to vary the relative weighting on the rigid body motion). In trans- 
forming the system to modal coordinates both K and M will be premultiplied 
T 

by T and postmultiplied by T, where T is the matrix of eigenvectors for the 
2nd order problem. Since the input influence coefficient matrix Bs’o (see 

T 

Table 4.2.1) is equal to T b, where b is a vector with a 1 in the first 
position and zeros in all other positions, the resulting Q matrix in 
modal coordinates will be: 


*2 = 


2 T 

A + q, B B 
n b o o 


] 


(4.2.57) 


For the purposes of this report the relative rigid body weighting factor q^ 

2 T 

has been set to 50. The resulting A + q^B o B Q matrix is listed in 
Table 4.2.3. 

Table 4.2.3 Potential Energy + Rigid Body Weighting Matrix 


.43 

-.57 

-1.33 

-2.38 

-.93 

-.52 

-.57 

49.07 

1.75 

3.13 

1.23 

.68 

-1.33 

1.75 

363.11 

7.28 

2.85 

1.59 

-2.38 

3.13 

7.28 

1350.40 

5.10 

2.84 

-.93 

1.23 

2.85 

5.10 

2694.60 

1.11 

-.52 

.68 

1.59 

2.84 

1.11 

3061.30, 
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Note that this matrix places considerably higher emphasis on the higher 

2 

modes. This is because the matrix A places a cost proportional to the square 
of the modal frequency. The other important difference between this weighting 
matrix and Q r is that a penalty is now placed on modal velocities. The closed 

loop pole configuration, as q c varies from 100 to 10,000, is illustrated in 

Fig. 4.7. An asterisk is placed next to the poles corresponding to 
q c = 1,000, and the loop gain for this case is plotted in Fig. 4.8. As 

expected, the control law resulting from places less emphasis on the rigid 

body modes, and more emphasis on the flexible modes, than does that resulting 
from Qj. This is indicated, both in the closed loop pole configuration, and 

in the loop gain. While LQR Design #1 pushed the rigid body poles further to 
the left than either of the first two sets of poles corresponding to flexible 
modes, LQR Design #2 pushes all the flexible mode poles further to the left 
than the rigid body poles. In comparing the two loop gains, that resulting 
from Design #1 shows a continual decrease in loop gain throughout the fre- 
quency range where flexible modes are being actively controlled, while the 
loop gain resulting from Design 2 indicates a non- decreasing loop gain 
through the same frequency range (though the gain is decreasing for both 
higher and lower frequencies). It is important to note though, that both 
designs actively control the flexible modes and have approximately the same 
bandwidth, though Design #1 does have a higher gain at low frequencies. Later 
in this Section it is found that the robustness properties of LTR designs 
based on these two different LQR designs vary considerably. The real rea- 
sons for this remain a topic for further investigation. 
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LQR Design #3 


There are two basic differences between Designs #1 and #2. The first is 
that Design # 2 places much more emphasis on the higher frequency modes, while 
the second is that Design # 2 places a penalty on modal velocities as well as 
modal displacements, while Design #1 does not. To separate these two effects. 
Design #3 uses the same penalties on modal position as #1 while adding identi- 
cal penalties to modal velocity. The resulting Q matrix is as follows: 




0 

Q. 


(4.2.58) 


The corresponding closed loop pole configuration as q varies from .01 to 10 

c 

is illustrated in Fig. 4.9. An asterix is placed next to the poles corres- 
ponding to q =1, and the loop gain for this case is plotted in Fig. 4.10. 
c 

Next results for asympototic LTR designs, corresponding to each of the LQR 
designs, will be presented. 


Asymptotic LTR Des igns 

There are a number of design variables in the asymptotic LTR procedure. 
These are the elements of the appended columns to the B-matrix, and a scalar 
parameter q, which is increased until loop recovery is achieved. Since only 
the case of two measurements is studied at this point there will be only one 
appended column to the B matrix. This is chosen as a scalar parameter r g 

times a vector corresponding to an input collocated with the second output 4 . 

4 The input influence coefficients corresponding to an input collocated with 
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The reason that a collocated input location is chosen is that this guarantees 

a minimum -phase- system for~x.-S-l.. The case r = 0 is equivalent to not 

s s 

adding any columns to the B-matrix, and does not strictly satisfy the condi- 
tions for asymptotic loop recovery outlined in Refs, [D01,2]. However, loop 
recovery is achieved for this case so it will be considered as a design possi- 
bility. 

Asymptotic LTR compensators corresponding to each of the three LQR designs 
are presented. In each case the loop gain for both the final step in loop 
recovery and the preceding step are plotted. 

The point at which loop recovery is achieved is somewhat of a relative mea- 
sure. As the parameter q is increased, three effects take place simultane- 
ously. First the loop shape begins to converge as pole/zero cancellations 
line up, second the overall loop gain increases and third the uncancelled com- 
pensator poles move out into the l.h.p. on 45° rays. Loop recovery is 
achieved when the overall shape and gain of the compensated plant loop match 
the full -state feedback loop and the decrease in loop gain due to the uncan- 
celled poles shifts past the point where loop gain has dropped below Odb. 

This will be clearer when the loop gain plots are presented in 

Figs. 4.11-4.17. The parameter q (listed as QEST on the figures) is actually 

multiplied by 10*^ to produce the results presented. This scaling was found 
to produce loop recovery for reasonable values of q, q is increased by steps 
of a factor of ten at a time, until it is judged, by visual inspection, that 


the 2nd measurement wil be identical to the output influence coefficients 
for the 2nd measurement, listed in Table 4.2.1. 


80 


loop recovery has been achieved. 


The tables describe robustness results. The first column gives the LQR 

design the second column the LQR parameter q , the third column the LTR 

c 

parameter q and the fourth column the 2nd input scaling factor r . The 

s 

robustness results indicate stability or instability of the closed loop system 
when all the frequencies are varied by the percentage indicated. An S indi- 
cates a stable closed loop system, while a U indicates an unstable closed loop 
system. This was chosen as a rough measure of robustness since the compensa- 
tors we tested seemed to be most sensitive to a shift of all the frequencies 
in the same direction. This is also a realistic possibility since a variation 
in mass or stiffness properties would produce just such a result. In a final 
analysis stage a more complex study of robustness (for all possible sets of 
frequency errors) would have to be carried out, but the data presented in 
Tables 4. 2. 5 -4. 2. 7 give an easily read and meaningful measure of robustness. 


Asymptotic LTR for LQR design #1 

In this case loop recovery was achieved for q=100. The loop gains for the 

case of q=10 and q=100 are plotted in Figs. 4.11 and 4.12 respectively. A 

number of different values of the scalar parameter r were tested and loop 

s 

recovery was achieved in every case. Robustness results are presented for two 

cases: r =1 and r =0. Neither case proved to be robust, but while the case 

s s 

with r =1 indicated stability only for a 10% increase in frequency, the case 
s 

with r -0 indicated stability only for a 10% decrease in frequency. This does 
s 
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suggest that the choice of the appended B-matrix columns will affect 
robustness. Note that increasing q to get loop recovery does not seem to 
improve robustness. In fact, the results do not change. However, increasing 
q does increase bandwidth and low frequency loop gain, thereby increasing per- 
formance while maintaing robustness. 

Table 4.2.4 Robustness of Asymptotic LTR Design for LQR Design #1 


Q# 

q c 

q 

r 

S 

-40% 

-30% 

-20% 

-10% +10% +20% +30% +40% 

l 

100 

10 

i 

U 

U 

U 

U 

S 

u 

u u 

l 

100 

100 

l 

U 

U 

U 

U 

S 

u 

u u 

l 

100 

10 

0 

U 

U 

U 

S 

U 

u 

u u 

l 

100 

100 

0 

u 

U 

u 

S 

U 

u 

u u 


Asymptotic LTR for LQR Design #2 

For the second LQR design, loop recovery was achieved for q=10. The loop 
gain for q=l and q=10 are plotted in Figs. 4.13 and 4.14 respectively, for 
this case different values of r g were not tested since the LTR desing was 

extremely robust for r =1. In fact Table 4.2.5 indicates stability for ±40% 

s 

variations in frequency at both steps in the loop recovery. This is an excel- 
lent robustness result, but we do not yet have a rigorous explanation for it. 
Table 4.2.5 Robustness of Asymptotic LTR Design for LQR Design #2 


0# 

% 

q 

r 

s 

-40% 

-30% 

-20% 

-10% +10% +20% +30% +40% 

2 

1000 

l 

i 

S 

S 

S 

S 

S 

s 

s s 

2 

1000 

10 

i 

S 

S 

s 

S 

S 

s 

s s 
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As indicated earlier, one difference between the two designs is that design 
#2 pushs all the closed loop poles corresponding to flexible modes further 
into the l.h.p. than poles corresponding to the rigid body mode, while design 
#1 does not. This indicates that the second design does not attempt to move 
the antenna rigid body mode too quickly, but it does pay a great deal of 
attention to the flexible modes. The first design, on the other hand, has a 
very high gain at low frequencies (indicating a fast response to errors in the 
rigid body mode), but does not pay much attention to the flexible modes. The 
performance of the second design, in terms of fast reduction in the RMS error 
or of the rigid body displacement error, will be less than that of the first. 
Though there is clearly a trade-off between performance and robustness, this 
does not entirely explain the difference in robustness, since much lower per- 
formance designs based on RMS error were also tried and no large increase in 
robustness was observed. This is still very much in the area of hand-waving, 
but a more rigorous understanding of the great increase in robustness observed 
for the second LQR loop shape, will be a major part of our effort in 1986. 


Asymptotic LTR for LQR Design #3 

LQR design #3 was concocted to separate the effect of the two major differ- 
ences between designs #1 and #2. That is the increased cost on the high fre- 
quency modal displacements and the cost on modal velocities. The modal dis- 
placement cost in design # 3 is identical to that of design #1, while a rate of 
change of RMS term is added to weight modal velocities. In this case, loop 
recovery was achieved for q=100. Loop gains for q=10 and q=100 are plotted in 
Figs. 4.15 and 4.16 respectively. Robustness results are presented in Table 
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4.2.6. This design is considerably more robust than #1, but not quite as 
robust as #2.— Again.,- we do not have a rigorous explanation of this fact, and 
this will remain an area of active research. Also note that increasing q does 
not im prove robustness, but it does improve performance. 

Table 4.2.6 Robustness of Asymptotic LTR Design for LQR Design #3 
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Algebraic LTR Designs 

The design variables in the Algebraic LTR approach are the locations of the 
compensator poles . The compensator order is fixed by the maximum degree of 
the plant numerator polynomials (n n ) and the number of measurements (m) . For 

the 6-mode antenna model with two measurements this will be order ten, or two 
orders less than an estimator based design. Of the ten poles available, eight 
will be cancelled and two will not. The positions of the uncancelled poles is 
set at -200±200j . This is far enough away from the controlled flexible fre- 
quencies that the loop shape in the design region is not significantly 
affected. To examine the effect of varying pole locations, four sets of pole 
locations were chosen. These are listed in Table 4.2.7. 
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Table 4.2.7 Compensator Pole Locations 


# 

Pole #1 

Pole n 

Pole #3 

Pole #4 

1 

-l±lj 

-212j 

-313j 

-414j 

2 

-lOHOj 

-20120j 

-30130j 

-40140j 

3 

-.15110 j 

-1.5125j 

-2. 5150 j 

-3. 0155 j 

4 

- 1.919. 8 j 

-5. 518. 3 j 

-8.315.5j 

-9.811.9 


The first two sets of poles lie on 45° rays from the origin. The next set 
are located close to the plant zeros in both channels. This is done in the 
spirit of the asymptotic approach which cancels plant transmission zeros. The 
fourth set of poles are placed in a Butterworth filter pattern with a cut-off 
frequency of 10 rad. /sec.. 

Loop recovery was satisfactorily achieved in every case, but robustness 
results were not particularily positive. Table 4.2.8 lists those results for 
all three LQR designs and all four sets of compensator poles. Only LQR design 
#2 in conjuction with pole location # 2 showed any respectable robustness 
results. Again, a rigorous interpretation of these results has not yet been 
made, but the number of possible cases for pole locations has not been 
exhausted. Once an ideal LQR Loop shape is found an optimizing routine might 
be applied to find the "best" pole locations. 
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Table 4.2.8 Algebraic LTR Designs 


Q# 

<0 

PL 

-40% 

-30% 

-20% 

-10% +10% +20% 

+30% 

+40% 

1 

100 

1 

U 

U 

u 

u 

u 

u 

u 

u 

1 

100 

2 

U 

U 

u 

u 

s 

u 

u 

u 

1 

100 

3 

U 

U 

u 

u 

u 

u 

u 

u 

1 

100 

4 

U 

u 

u 

u 

u 

u 

u 

u 

2 

1000 

1 

U 

u 

u 

u 

u 

u 

u 

u 

2 

1000 

2 

S 

s 

s 

s 

s 

s 

u 

u 

2 

1000 

3 

u 

u 

s 

s 

u 

u 

u 

u 

2 

1000 

4 

u 

u 

u 

u 

u 

u 

u 

u 

3 

1 

1 

u 

u 

u 

u 

u 

u 

u 

u 

3 

1 

2 

u 

u 

u 

u 

u 

u 

u 

u 

3 

1 

3 

u 

u 

u 

u 

u 

u 

u 

u 

3 

1 

4 

U 

u 

u 

u 

u 

u 

u 

u 


92 


4.3 SENSITIVITY OPTIMIZATION 


Our second approach to robust compensator design is to use parameter optim- 
ization (nonlinear programming) to minimize various measures of the sensitiv- 
ity of the closed-loop system to parameter errors in the model. Initially, we 
augmented the quadratic performance index in the optimal linear regulator 
problem with terms penalyzing the gradient of the closed-loop response (sensi- 
tivity) with respect to uncertain parameters. While the resulting optimal 
control problem could not be solved via the solution of a Riccati equation, 
the function and gradient evaluations for parameter optimization involved 
solution of Liapunov matrix equations, which can be solved rather efficiently 
by well known methods. Unfortunately, this approach proved ineffective 
because we were unable to find a quadratic performance index involving sensi- 
tivity vectors such that the compensator that minimized the performance index 
was robust. 


A much more successful approach has been to minimize an objective func- 
tional involving the gradients of the real parts of the closed- loop eigenva- 
lues with respect to uncertain parameters; i.e., minimize the sensitivity of 
the closed-loop eigenvalues. The closed-loop system has the form 


fl 



(z) 


<8 


where 
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-BK 

(A-BK-GC) 

Our typical sensitivity objective functional has the form 
J s = Z[dRe(X i )/da j ] 2 

where the X^*s are the eigenvalues of A^j and the a^'s are the uncertain par- 
ameters in A, B and C. The uncertain parameters usually are natural frequen- 
cies and damping ratios. The design variables for this optimization problem 
are the control and estimator gains K and G. 

Usually, we obtain the initial gains as the solution to an optimal LQG 
problem, and then use sensitivity optimization to improve the robustness of 
the closed-loop system. We use standard formulas [KA2,Ch. II .2] for computing 
dX^/da^ in terms of the eigenvectors and eigenvalues of A^ and dA^/da^, and 

we have found that the computation is much more efficient if done entirely in 
real arithmetic. 

For our beam-hub model, we have been able to take an initial LQG compensa- 
tor that allows less than 10% variation in the natural frequencies before the 
closed-loop system becomes unstable and obtain a sensitivity-optimized compen- 
sator that allows more than 20% variation in the frequencies while maintaining 
closed-loop stability. 

The following two tables list the dominant closed- loop eigenvalues of the 
antenna and two compensators based on a six- mode antenna model. The first 
compensator is optimal for an LQG problem in which state weighting Q penalizes 
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the r.m.s. surface error of the antenna and the disturbance is white noise 


distributed uniformly over the antenna, along with white noise entering 
through the actuator. This compensator allows less than 10% error in the 
open-loop plant frequencies before the closed- loop system becomes unstable. 
With the LQG compensator as the initial guess, our sensitivity optimization 
program found the second compensator, which allows more than 20% error in the 
open-loop frequencies while maintaining closed-loop stability. 


95 



Table 4 .3.1. Dominant Closed- loop Eigenvalues 
with LQG Compensator 


For correct frequencies For +10%' error in 

all frequencies 
-1.2 ± ,2j -2.0 


For +20% error 


-10.7 


- 1.1 


1.0 


-3.3 ± 4. 8j 
-3.0 ± 6 . 7 j 
-1.1 ± 7 . Oj 
- .6 ± 19. j 
-2.8 + 20. j 
-3.8 + 30. j 


-4.8 ± ,5j 

- .4 ± 6 . Oj 
-2.7 ± 7 . Oj 

+ -4 ± 17. j 

- 3.5 ± 20. j 

- .5 ± 33. j 


- 1.6 ± 1.3j 

- .06 ± 5.6j 
-2.50 ± 7 . Oj 

+ .80 ± 16. j 

- 3.8 ± 19. j 
+ 0.2 ± 36. j 


Table 4.3.2. Dominant Closed-loop Eigenvalues with 
Sensitivity Optimized Compensator 


For correct frequencies 
-1.3 


For +10% error in 
all frequencies 
- 1.6 


For +20% error 
- 2.2 


.6 


.6 


.6 


-3.4 ± 4.8j 

- .6 ± 6 . 9 j 
-3.0 ± 6 . 7 j 
-2.8 ± 19. j 

- .7 ± 20. j 
-2.2 ± 36. j 


-4.1 ± 4. Oj 

- .6 ± 6.0j 
-2.5 ± 6 . 7j 

- .3 ± 17. j 
-3.2 ± 20. j 
-1.3 ± 33. j 


-4.5 ± 2.5 j 

- .6 ± 5 . 2 j 
-2.1 ± 6. 8j 
-.05 ± 16. j 
-3.6 ± 20. j 

- .7 ± 29. j 
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We have found that the form of the operator Q in the optimal control prob- 
lem affects robustness significantly. The Q for the preceding LQG compensator 
penalized the rigid-body displacement primarily and did not penalize energy at 
all. Table 4.3.3 gives the dominant closed-loop eigenvalues produced by the 
LQG compensator for a Q which penalizes both the rigid-body angle and the 
total energy in the antenna. We chose the Q-to-R ratio so that the real part 
(-.65) of the least stable closed-loop eigenvalue with this compensator is 
comparable to the real part (-.6) of the least stable eigenvalue with the pre- 
vious LQG compensator (which resulted from the Q that penalized the antenna 
surface error). In this case, sensitivity optimization did not improve the 
robustness of the LQG compensator significantly. 

Table 4.3.3. Dominant Closed- loop Eigenvalues 
with Second LQG Compensator 

Ed correct frequencies For +10% error in For +20% error 


all frequencies 


-0.65 

+ 

0. 13j 

-0.56 

± 

0. 17 j 

-0.53 

+ 

0.18 j 

-0.70 

± 

0.33j 

-0.76 

± 

0.42 j 

-0.76 

+ 

0.48 j 

-0.78 

± 

6.94j 

-1.36 

± 

6.83j 

-1.50 

+ 

6.92 j 

-0.84 

± 

6.96j 

-0.27 

± 

7.81 j 

-0.13 

± 

8.46 j 

-1.13 

± 

18. 9 j 

-2.49 

± 

18. 8j 

-2.61 

± 

18.9 j 

-2.05 

± 

18. 9j 

-0.71 

± 

21. lj 

-0.59 

± 

22.8 j 
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4.4 OPTIMIZED LOOP RECOVERY DESIGNS 


Another approach to improving robustness combines the Loop Transfer Recov- 
ery methods described in Section 4.2 with the sensitivity optimization methods 
described in Section 4.3. This is the optimal Loop Transfer Recovery 
approach, and is currently under investigation. 

The idea of the approach is to take advantage of the robustness improvement 
inherent in the sensitivity optimization methods, while maintaining a fixed 
loop shape. This guarantees that robustness is not achieved at the cost of a 
reduced loop gain, or equivalently reduced performance. As noted in Section 
4.2.2, the loop recovery compensator for a given full state feedback loop 
shape is not unique when a plant has more outputs than inputs . In the asymp- 
totic LTR approach, this non-uniquenss corresponds to the choice of artificial 
columns in the B-matrix. This allows the designer some freedom in the place- 
ment of finite plant transmission zeros, which in turn determine the location 
of the compensator poles. In the algebraic LTR approach the degrees of free- 
dom depend on a number of factors including the desired compensator order, 1 
but for the specific algorithm described in Section 4.2.3, the compensator 
order is fixed and the degrees of freedom translate directly into the location 
of the compensator poles. 

Consequently two optimal LTR approaches are suggested. The first finds an 
optimal set of elements for the artificial columns of the B-matrix in an 
asymptotic LTR design. The second finds an optimal set of compensator pole 


1 As compensator order is increased, the number of degrees of freedom in the 
algebraic loop recovery design also increases. 
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locations in an algebraic LTR design. In both cases the sensitivity of the 
real parts of each closed loop pole to variations in plant frequency are eval- 
uated. The objective function is equal to the sum of the squares of these 
sensitivities. This objective function is used since it gave excellent 
results in the more general sensitivity optimization approach described in 
Section 4.3. 

Since the data presented in Section 4.2.5 suggests that the asymptotic LTR 
designs provide greater robustness than the algebraic designs, the optimal 
approach is being applied first to the asymptotic design method. Currently 
the software is written and preliminary results are expected early in 1986. 

Two other applications of optimization to the LTR approach are also being 
considered. The first is based upon the conjecture that the lack of robust- 
ness exhibited by LTR designs relative to the equivalent LQR designs, is a 
function of pole/zero cancellations close to the imaginary axis. As noted in 
Section 4.2.4, small variations in plant parameters can result in Very large 
variations in loop shape whenever pole/zero cancellations occur close to the 
imaginary axis. An optimization routine could therefore be used to choose 
artificial columns of the B-matrix so as to shift finite plant transmission 
zeros away from the imaginary axis. A second possibility depends on the fact 
that the LTR estimator design is largely independent of the full-state feed- 
back LQR control gain matrix K. An optimization routine could therefore leave 
the LTR estimator design fixed while varying the elements of K so as to minim- 
ize closed loop pole sensitivity. Unlike the other optimal approaches 
described in this Section, this approach would not maintain a given loop 
shape, and performance might be reduced, but it may produce some further 
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insight into the relationship betweeen full-state feedback loop shape and LTR 
robustness . 

These optimal LTR approaches to achieving robust control of the antenna 
will be studied more extensively in 1986. 
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5.0 CONTROL DRIVEN FINITE ELEMENTS 


5.1 A COMPARISON OF THREE FINITE ELEMENT SCHEMES FOR 
APPROXIMATING THE SOLUTION TO THE OPTIMAL LINEAR 
REGULATOR PROBLEM FOR A FLEXIBLE STRUCTURE 

The Control System 

One end of the Euler-Bernoulli beam in Figure 5.1.1 is attached rigidly 
(cantilevered) to a rigid disc which is free to rotate about its center, point 
0, which is fixed. Also, a point mass m^ is attached to the other end of the 

beam. The control is a torque u applied to the disc, and all motion is in the 
plane. 
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Table 5.1.1. Structural Data 


r = hub radius 
1 = beam length 

I Q = hub moment of inertia about axis 

perpencilar to page through 0 
m^ = beam mass per unit lingth 


10 in 
100 in 

100 slug in 2 
.01 slug/ in 


m^ = tip mass 1 slug 

El = product of elastic modulus and 

second moment of cross Section for beam 13,333 slg in 3 /sec 2 

fundamental frequency of undamped structure .9672 rad/sec 


The angle 0 represents the rotation of the disc (the rigid- body mode) , 
w(t,s) is the elastic deflection of the beam from the rigid-body position, and 
Wj(t) is the displacement of m^ from the rigid-body position. The control 

problem is to stabilize rigid-body motions and linear (small) transverse elas- 
tic vibrations about the state 0=0 and w = 0. Our linear model assumes not 
only that the elastic deflection of the beam is linear but also that the axial 
inertial force produced by the rigid-body angular velocity has negligible 
effect on the bending stiffness of the beam. The rigid-body angle need not be 
small. 

This system can be written in the form (3.1.1) with the generalized dis- 
placement vector 

x = (0,w,w 1 ) H = R X L 2 (0,1) X- R. (5.1.1) 
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We assume that the beam has Voigt-Kelvin viscoelastic damping, so that the 
damping operator in (3-1.1) is 


D 0 ~ C 0 A 0 


(5.1.2) 


where c^ = 10 . For the details of the energy-space formulation of this 

problem and the operators in (3.1.1), see [611]. 

The Optimal Control Problem 

In the optimal control problem, we take Q = I in the performance index in 
(3.2.1). This means that the state weighting term <Qz,z>g is twice the total 

energy in the structure plus the square of the rigid-body rotation. Since 
there is one input, the control weighting R is a scalar, and we take R = .05. 
The functional control gains in (3.3.2) have the form 

f = (a f ,0 f ,0 f ), (5.1.3) 

S = ( , PjT ) J (5.1.4) 

where the a's and 0's are scalar gains for the hub and tip mass and 0. and 0 

f g 

are functions which define the feedback law for the beam. Because the strain 
energy in the beam involves the second spacial derivative of w, 0^ t , appears 

in the the control law in (3.3.2), rather than 0^. 

Approximat ion 


To approximate the solution to the optimal control problem, we solve a 
sequence of finite dimensional optimal control problems as described in Chap- 



ter 3.4. To illustrate the convergence of the approximating feedback control 

laws, we will plot jthe_dis.tribu£ed functional gain components $J' and $ . 

f 8 

Ve compare three approximation schemes; i.e., three different sets of basis 
vectors e^ for the approximation to the generalized displacement in (3.4.1). 

The first approximation of the structure is based on a finite element approxi- 
mation of the beam which uses cubic Hermite splines [ST1] as basis functions, 
the second approximation scheme uses cubic B-splines to approximate the beam, 
and the third approximation scheme uses the natural mode shapes of the struc- 
ture as the basis vectors in the approximating optimal control problems. 

Recall that cubic B- splines require the displacement, slope and second deri- 
vative to be continuous at the nodes, while cubic Hermite splines require only 
the displacement and slope to be continuous at the nodes. To get the natural 
modes of the structure, we used 24 elements in the Hermite spline approxima- 
tion to get the first 12 modes of undamped, free vibration. 

Figure 5.1.2 shows the distributed functional gain components computed 
using the Hermite splines. Figure 5.1.3 shows the distributed functional gain 
components computed using the B-splines, and Figure 5.1.4 shows the distri- 
buted functional gain components computed using the natural mode shapes. For 
both spline schemes, we get convergence with ten elements and very near con- 
vergence with eight elements. Also, we get convergence with ten modes and 
very near convergence with eight modes. 

Now note that if n @ is the number of elements, then the number of degrees 
of freedom in the Hermite spline scheme is 2n @ + 1, while the number of 
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degrees of freedom in the B-spline scheme is n + 1. This means that the 

B-splines and natural mode shapes are almost equally efficient for computing 
the control gains and that the Hermite splines require almost twice as many 
degrees of freedom. This is important because the dimension of the Riccati 
matrix equation that must be solved in each approximate optimal control prob- 
lem is twice the number of degrees of freedom in the approximation scheme. 

We experimented with varying the location of the nodes for the Hermite 
splines and B-splines. While we found that we could improve the rate of con- 
vergence slightly, we felt that the improvement was not sufficient to justify 
the elaborate procedure. The improvement in convergence rate obtained by var- 
ying the node locations was insignificant compared to the improvement obtained 
by using modes or B-splines instead of Hermite splines. 
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Figure 5.1.2a 
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Figure 5.1.2b 
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Figure 5.1.3b 
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Figure 5,1.4b 
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5.2 DIRECT APPROXIMATION OF HAMILTONIAN SYSTEM 


We have investigated approximation of the solution to the optimal control 
problem for the distributed model of the space stucture by constructing direct 
finite element approximations of the infinite dimensional Hamiltonian system. 
To date, our numerical solution of the LQG optimal control problem for the 
distributed model of the space antenna has been based on an approximation 
scheme that is indirect in that we first approximate the open- loop system with 
a sequence of finite dimensional systems and then solve a finite dimensional 
optimal control problem for each approximating system. The solutions to these 
finite dimensional problems converge to the solution to the infinite dimen- 
sional problem. However, for certain space structures, direct approximation 
of the solution to the infinite dimensional problem may allow better choices 
for the basis elements used in the approximation. 

The Hamiltonian operator involved in the solution to the optimal regulator 
problem is 



where A is the open- loop semigroup generator, B the input operator, and Q the 
state-weighting operator in the optimal control problem. Recently, we began 
studying the extension of an idea used at JPL for solving static state estima- 
tion problems. As applied to the control problem, the idea is to choose a 
finite element scheme to approximate the infinite dimensional Hamiltonian 
operator directly with operators of the form 


PRECEDING PAGE BLANK NOT HMD 
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(5.2.2) 


X 

Here, is not necessarily the adjoint of A^. In particular, two different 

* 

sets of basis vectors may be used in approximating A and A . 


For the models of flexible structures and the finite element schemes that 
we have used to date, direct approximation yields the same finite dimensional 
Hamiltonian operator as does indirect approximation, but for more complicated 
structure models, this should not be the case. Direct approximation of the 
infinite dimensional Hamiltonian operator may be most useful for structures in 
which some of the materials are viscoelastic with history-dependent damping. 
For discussion of history-dependent viscoelastic effects see [C01,FU1]. Such 
structures lack the basic selfadjointness of more common structures, and the 
* 

operators A and A have very different domains. Hence different finite ele- 

* 

ment basis vectors are needed for approximating A and A . These basis vectors 
should be chosen to produce appropriate convergence of the direct approxima- 
tions to the Hamiltonian system and the corresponding direct approximations to 
the solution of the infinite dimensional optimal control problem. 
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6.0 SUMMARY AND CONCLUSIONS 


The work described in this report further develops aspects of modern con- 
trol theory which make it better able to deal with issues of modeling and 
implementation. The report represents an extension of previous work done for 
JPL in which elements of distributed parameter control theory and finite 
dimensional control theory have been synthesized into a design approach which 
simultaneously produces both a model which is appropriate for the control 
problem at hand, and an "optimal" controller based on that model. The inter- 
action of modeling and controller design is central to the approach. The 
techniques developed have been applied to a space antenna to demonstrate their 
performance on space systems of realistic complexity. 

A major focus of the work carried out during 1985 was robustness. Robust- 
ness refers to the ability of a control system to perform satisfactorily even 
when the model used for the design is an imperfect representation of the 
actual system. There are two types of modeling errors which are of major con- 
cern in designing control systems for large space structures (LSS). Much of 
the work on functional gains is aimed at one particular type of modeling prob- 
lem: model truncation. Parameter errors are another important source of 
modeling error. Work during the current year emphasized the development of 
LSS controller designs which are relatively insensitive to parameter errors 
but still maintain a satisfactory level of system performance. For antenna 
structures of the type considered, this is a difficult problem which does not 
easily yield to standard techniques. 

As a general rule, an improvement in robustness can only be achieved at the 
cost of a reduction in some measure of performance. In order to achieve a 
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proper balance between performance and robustness, some constraints must be 
imposed. In the loop recovery approach, the open loop and closed loop fre- 
quency responses are constrained to approximate those of the full state feed- 
back LQR design. Thus, the performance of the system is first determined in 
the context of the LQR design, and then one tries to obtain an estimated state 
feedback implementation which has similar performance and also adequate 
robustness. The sensitivity optimization approach also begins with the solu- 
tion of an LQR problem and designs an estimated state feedback implementation. 
In this case, however, there is no constraint to retain the LQR frequency res- 
ponse. Instead, the goal is to minimize the sensitivity of the closed loop 
eigenvalues to parameter errors. A constraint must be imposed to keep the 
closed loop regulator and estimator eigenvalues sufficiently far in the left 
half plane otherwise the optimization routine would generate a very low per- 
formance control system. Designs based on sensitivity optimization tend to be 
more robust than LTR designs, but their loop gains are also lower. 

Both of these approaches rely heavily on the statement of the LQR problem 
which serves as a starting point for the designs. The state weighting opera- 
tor in the performance index for the LQR problem discussed in Chapter 3 is Q. 
Generally, scaling up Q while holding the control weighting matrix R fixed 
leads to larger control gains and a less robust compensator. The numerical 
examples presented in Sections 4.2 and 4.3 indicate that in addition to the 
magnitude of Q, the relative penalization that Q specifies among the different 
components of the state vector can affect robustness dramatically. For exam- 
ple, Q's which penalize only the RMS surface error of an antenna model were 
found to lead to designs which were not very robust. However, Q's which 
penalize the energy or the RMS surface error plus the rate of change of RMS 
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error, were found to be quite robust. While the interpretation of these 
results is not entirely clear at this point, some observations can be made. 
Penalizing RMS surface error alone leads to controllers which try to reduce 
rigid body displacement much faster than elastic deformations. This means 
that the closed loop eigenvalues for the rigid body mode are substantially to 
the left of those for the lowest flexible modes. In this case, large control 
torques are produced to control the rigid body mode, and the system is not 
robust because these torques strongly excite the flexible modes which may be 
inaccurately known. The other two Q' s mentioned above penalize the rigid body 
and flexible modes more equally. The real parts of of the closed loop eigen- 
values, and hence the decay rates, for the rigid body mode and lowest elastic 
modes are approximately equal and the compensator is more robust. In our con- 
tinuing research effort, we plan to pursue further the relationship between 
the selection of Q and robustness. 

Another factor which has a strong effect on robustness in the number an 
position of the measurements. In the case of the beam-hub model, control 
designs based on two measurements were much more robust than those based on 
one measurement. The implications of the number of measurements for the 
antenna model will be studied more extensively in 1986. 

Sensitivity optimization appears to be a successful approach for designing 
robust compensators for complex structures, but we need more analysis and num- 
erical experimentation to determine the best objective functionals and the 
most effective methods of optimization. We have experimented with objective 
functionals involving both first and second-order gradients of eigenvalues 
with respect to uncertain parameters, and found that minimizing the second 
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gradients also does not improve robustness significantly. Also, we are 
comparing performance of the optimization routines using analytic as opposed 
to numerical gradients with respect to the design variables. 

The results in Section 5.1 demonstrate that a well chosen finite element 
scheme can greatly reduce the order of approximation required to design an 
optimal compensator for a flexible structure. It should not be surprising 
that the natural mode shapes are at least among the best basis vectors for 
approximating the optimal functional gains. It is more interesting that the 
cubic B -splines are almost as efficient in the example that we worked; they 
are much better than the cubic Hermite splines, which are most common in engi- 
neering finite element approximations of beams. This suggests that further 
research on control driven finite elements could produce significantly more 
efficient numerical methods for the design of controllers for large flexible 
structures . 

Since we know now that some approximation schemes are more efficient than 
others for numerical solution of distributed optimal control problems, it 
seems natural to try to develop a scheme especially for the approximate solu- 
tion of the optimality conditions for the infinite dimensional LQG problem for 
the distributed model of a flexible structure. These optimality conditions 
are embodied in the infinite dimensional Hamiltonian operator discussed in 
Section 5.2, where we propose the direct approximation of this operator 
instead of the indirect approximation that usually is used. For the kinds of 
material models and finite element schemes that we have used so far, the 
direct and indirect approximations are identical, but for materials with 
memory (particularly history-dependent damping) the direct approximation of 
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the Hamiltonian operator would be different and should be investigated. 
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Appendix A: Singular Values 


Singular values measure the "size" of a complex matrix, and are equivalent 
to the concept of the modulus of a complex scalar. The singular values of a 
matrix transfer function can therefore be thought of as the MIMO equivalents 
of the gain of a scalar transfer function. 

Any n x m complex matrix A can be decomposed in the following fashion: 


A ^ 

A = UIV = X o.u.v. 

£l 111 


* 


where : 

k = min(n,m) 

o. =A7 Taa ]=A7[a a] , 


°i 4 °i + l 



U = [.J.J ... u n ], 


* 

uu 

= U U = I 

V = [Vj v 2 ... vj. 

* 

* 

w 

= V V = I 


The o/s are the singular values of A, the vectors u^ are the left singular 
vectors of A and the vectors v^ are the right singular vectors of A. The 
largest singular value is denoted by 7 and the smallest singular value 
by o. 7[A] is the 1^ induced operator norm of A in the following sense: 
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TJ[A] = max ||Ax||, 


x^O 


lx « 


max || Ax || _ 

II x » 2 =1 2 


while o_[ A] is a measure of the closeness of A to singularity in the following 
sense: 


£[A] = min ||Ax|L = 

xs*0 I, H 
llx|L 


min || Ax || 

llx 1 2 =1 2 


&[A] is the only computationally reliable tool for the determination of near 
singularity, or rank of a matrix [KL1,LA1]. 


The following table includes a number of useful properties of singular 
values : 


1) £[A] > 0 a[A] = 1/7JIA" 1 ] 


2) 

^[“A] =| oe 1 0 i [A] 



3) 

o[A] 5 |X.[A]| < 

o[A] 


4) 

A=A* o i [A] = 

|X i [A]| 


5) 

A=A*£0 T~~^ o i [A] 

= X.[A] 


6) 

o[A + B] S 3[A] - 



7) 

■3[AB] £ "o[A]o‘[BJ 



8) 

JltAB] £ £[A]a[B] 



9) 

SL[A] - 1 5 o[I + 

A] £ oJA] 

+1 

10) 

V[E] < o[A] < 

J[A + E] > 

0 

11) 

TT[A] < 1 I 

+ A] £ 1 ■ 

-TT[A] 
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Another slight variation on the singular value decompostion is the polar 
decompostion, used extensively by Postlethwaite [P01,2,K01]. Write the singu- 
lar value decompostion as: 

A = UIV* = UV*(VXV*) = (UZV*)UV* 

it ie it 

where UZV is the usual singular value decomposition, UV (VZV ) is the right 

it it 

polar decomposition and (UIU )UV is the left polar decomposition, Postleth" 

it it 

waite defines the principal gains as the eigenvalues of VZV or UZU * which 
are of course just the singular values of A, and the principal phases as the 

it 

arguments of the eigenvalues of the unitary matrix UV . This decomposition is 
useful in separating a transfer function into a gain part and a phase part and 
will be used the proof of Theorems 2.1 and 2.2 in Appendix B. 
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Appendix B: Proof of Unstructured Uncertainty Theorems 

Theorem 4.1 will be proved in detail. The changes that need to be made for 
Theorem 4.2 will not be extensive and will be described in less detail. 

Theorem 4.^: 

The Multivariable Nyquist Criterion requires that det[I + KG'], evaluated 
on the standard Nyquist D-contour, encircle the origin counterclockwise as 
many times as KG' has unstable open-loop poles. If it is assumed that KG' has 
the same number of unstable open- loop poles as KG, the encirclement count of 
the origin must not change as G is warped continuously towards any allowable 
G' . For the Nyquist Contour to remain fixed as G is warped towards G' , it is 
also necessary to assume that any poles on the jw-axis be identical for G and 
G 1 . Any plant that has an uncertain frequency on the jw-axis will have an 
infinite unstructured error in the neighborhood of that frequency, so it is 
reasonable to assume that the only poles on the jw-axis are at the origin, or 
equivalently, that the model have some positive damping. Furthermore, practi- 
cal plants will be strictly proper, implying that lim K(s)G(s)=0. The Nyquist 

S"*’«» 

D-contour can therefore be replaced by the jw-axis, with the possibility of an 
indentation about the origin, which will not have any effect on the results 
for uncertain modal data. 

Since KGCI+eA^) is a continuous function of c, requiring that the encircle- 
ment count of det[I + KG] not change as G is warped towards G' is equivalent 
to requiring that: 

detfl + KG(I+eA )] * 0 for 0 £ e £ 1 
m 
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or in terms of the minimum singular value: 


j 


o,[I + K(s)G(s) (I+eA (s))] > 0, for 0<s£l, o-[A(s)] < 1 (s), seA., (B.l) 

m in K 

Dropping the dependence on s, and noting that if , o[A(s) ]-l m > then TT[eA(s) ]- l m > 
allows us to replace (B.l) by: 

£[I + KG +KGAJ > 0 for 7[A] £ 1 (B.2) 

Assume that fi.[(I + KG) ^]=1/C7[I + KG]) > 0. This is true whenever 
7[KG] < •», which in turn is true by the definition of the Nyquist D-contour. 

Multiply (B.2) byjd[(I + KG)" 1 ]: 

0 < o[(I+KG) _ 1 ]2.[I + KG +KGA ] < jl[I + (I + KG) _1 KGA ] (B.3) 

m m 


Now let (I 4* KG) *KG have the singular value decompost ion given by UZV or 
equivalently the polar decompost ion given by UV (VIV ). Let A = -el^VU for 
O^e^l. ii[ A] = 1T[A] = 1 1^ £ 1^, so this is an uncertainty that falls within 
the specified bounds. Substitute this in: 

o < a[I + (I + KG)" 1 KGA ] * J2[I - El VIV*] (B.4) 

m m 


but (I-eI^VXV ) is hermitian, positive semi -definite, so the singular values 
can be treated as eigenvalues to get: 


0 < 1 - El m o i (VXV*) = 1 - el m O i [(I + KG) _1 KG] 


(B.5) 
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(B.6) 


1 - *1 V[(I + KG) -1 KG] > 0 

*oT[ (I + KG) -1 KG] > 1/ Cel ) 1/(1 ) (B.7) 

mm 

Before continuing to prove sufficiency it is worth noting that the particu- 
lar uncertainty that will destabilize the system must have a phase character- 

ie 

istic given by -VU . An uncertaintly with some other phase characteristic 

might require that o[A (jw)]»l (w) before the system would go unstable. This 

mm 

indicates that if there is a some specified relationship between gain uncer- 
tainty and phase uncertainty, the unstructured uncertainty representation of 
Theorem 1 can be arbitrarily conservative. This point is discussed by Lehto- 
maki [LEI, 3, 4] in the context of most sensitive directions of the perturba- 
tion. 

Now prove sufficiency by working backwards: 

~[(I + KG) _1 KG] < 1/(1 ) (B . 8) 

m 

ZT^l > l m o[(I + KG) _1 KG] £ TT[A]T[(I + KG) _1 KG] (B.9) 

+ (I + KG)" 1 KGA m ] S 1 - “Old + KG) _1 KGAJ > 0 (B.10) 

Now assume that £[I + KG] > 0, which is true whenever the nominal feedback 
system is stable. Multiply by £[1 + KG]: 

0 <*[! + (I + KG) _1 KGA m ]i[I + KG] ^i[I + KG + KGAJ (B.ll) 
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(B. 12) 


o[I + KG + KGA ] > 0 
m J 

and Theorem 4.1 is proved 1 . 

Theorem 4.2: 

The major difference in the proof of Theorem 4.2 is in ascertaining the 
continuity of L(e,s) as s varies from 0 to 1. In particular KGCI+cA^) * must 

be a continuous function of e for 0£eSll, or equivalently L(s)=(I+A^) * must 

have no zero or strictly negative real eigenvalues. This is true since if 

L(s) has no zero or negative eigenvalues, then neither does (I+A.), and thus 

a 

A d can have no eigenvalues in the interval (-“,-1], so eA^ never has eigenva- 
lues of -1 and (I+eA^) is never singular. If the first requirement is met it 
is then necessary that: 

0 <a.[I + KG(I+A d )" 1 ] Vseft R , Vee[0,l] (B.13) 

which in turn is equivalent to: 

0<jj[I + KG + A d ] Vseft R (B.14) 

The proof then continues as with Theorem 4.1. 


1 The only difference between our result and the result stated in Ref. [D01] 

is that we've replaced ‘of A(jw)]<l by 7[A(ju»)]Sl . This is a minor varia- 

no m 

tion, and it was corrected in later papers by Doyle [D03,4,5]. Furthermore, 
for practical purposes, the results are identical. 
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Appendix C: Numerical Considerations in Algebraic LTR 

There are essentially five steps in solving for an LTR compensator by the 
algebraic method. This appendix will describe the numerical considerations 
important in each step. 

1) Calculation of the Plant Numerator Polynomials - 

The plant transfer function is as follows: 

G(s) = C(sI-A)” 1 B (C.l) 

where C is an mxn matrix, A is an nxn matrix, B is an nxl vector and G(s) is 
an mxl vector function of s. To use the algebraic approach the numerator 
polynomials of G(s) must be found. We have used two methods, the Fadeeva 
method, which finds the coefficients directly based upon the Cayley Hamilton 
Theorem [CHI], and a method based on the transfer function zeros. The Fadeeva 
method works well for small systems but not for large. The first reason for 
this is that the number of steps is proportional to n*, which can become com- 
putationally very expensive. The second reason is that the algorithm involves 
raising the matrix A to the (n-l)th power. This works well for discrete time 
systems, since the eigenvalues of A lie inside the unit circle, but it is very 
poorly conditioned for continuous time problems, since the eigenvalues of A 
may be very large. 

The second approach is far more reliable. It is based on calculating the 
zeros in each channel of the plant, calculating a constant that will premulti- 
ply the numerator polynomial and then multiplying out the zeros. The zeros of 
a plant are defined as values of s which result in a zero output for a non- 
zero input. In terms of the state space representation of the plant for the 
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ith channel, this can be defined as follows: 




(sI-A) -B 

H -W 

C. 0 

i J 

M b 


where is the ith row of C and the plant zeros are the values of s which are 

a solution to (C.2) for y = 0. This can in turn be transformed into the fol- 
lowing generalized eigenvalue problem: 


r a b 

1 

H 

= S 

I 0 

1 

[j 

i 

n 

o 

J 

fu| 

l 

0 0 

J 

u 

i j 

< 

k 4 

) 


< 



(C.3) 


The generalized eigenvalue problem can be solved accurately and efficiently by 
the QZ - algorithm [ST2] , where the computational cost is proportional to n k . 

The only tricky step is the separation of finite zeros from infinite zeros. 

The QZ-algorithm will result in one complex vector a and one real vector fJ. 

The eigenvalues are found by dividing the elements of a by the corresponding 
non-zero elements of f$. The first step in separating finite zeros is to eli- 
minate all elements where f$ is zero. There are two possibilities for elimi- 
nating further infinite zeros. The first is to place a threshold on the ele- 
ments of (S, eliminating elements lying below this threshold. This method, 
however, proves to be very sensitive to the threshold. A more reliable method 
involves dividing all the remaining elements of o by the corresponding non- 
zero elements of and then placing an upper bound on the magnitude of the 
resulting zeros. The upper bound is easily estimated, with some knowledge of 
the problem, and this method has proved to work very well. As a final check, 
the number of zeros for a SISO flexible system will almost always be n-2. 


The ith numerator polynomial will have the following form: 


131 



t 


n^(s) = k^s-z^s- z 2 ) • . • (s-z n _ 2 ) 


(C.4) 


The next step therefore involves computing k^ This is most easily done by not- 
ing the following equality: 


, k (s-z )(s-z ) . . . (s-z .) 

C.(sl-A) B = ~ 77 — — " 

x (s-p 1 )(s-p 2 )...(s-p )(s-p )(s-p ) 


(C.5) 


where Pj is the jth pole. (C.5) is valid for any s which is not a pole of the 
system, k^ can therefore be found by evaluating (C.5) for some s = s q 

(eg. s q = 1) and then solving for k^. Also note that only (s q I-A) *B need be 

evaluated, using a linear equation solver, rather than (s I -A) *. The above 

o 

method was used in Ref. [EDI] to evaluate plant transfer functions for s = jw. 

The final step in solving for the plant numerator polynomials is to multi- 
ply out (C.4). The approach based on plant zeros is more complex than the 
Faddeva method and must be carried out for each channel separately, but it is 
far more computationally reliable. 

2) Solution to the LQR Problem - 


The major step xn the calculatxon of full— state feedback control gaxns is 
the solution of a Riccati equation for the LQR problem. This is solved by 
Potter's method [KA1] which involves finding the eigenvalues and .eigenvectors 
of the following Hamiltonian matrix: 


H = 


-1 T 
-BR V 


-Q 


-A 


(C. 6 ) 


Finding an accurate solution to the Riccati equation therefore reduces to 
finding accurate eigenvalues and eigenvectors for (C. 6 ). Inaccuracies can 
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develop when Q is made very large in comparison with BB . This becomes espe- 
cially important in the solution of the near singular Kalman filter problem 
for the asymptotic LTR method, since Q may be very large before loop recovery 

T 

is achieved. This difficulty can be overcome by multiplying both Q and BB by 
the square root of q^, and then multiplying the final Riccati solution matrix 

P by the same number, rather than simply multiplying Q by q^. This procedure 

keeps H from becoming too "unbalanced" and results in accurate solutions to 

the Riccati equation, even for very large q . 

c 

The numerator polynomial for the optimal LQR transfer function K(sl-A) *B 
is found by the same methods used to find the plant numerator polynomials. 

3) Compensator Order Determination - 

The compensator order is found very simply by counting the maximum degree 

of the plant numerator polynomials and applying equation (4.2.24). If a 

proper compensator is required the compensator order will be identical to n VT 

N 

in (4.2.24), while if a strictly proper compensator is required, the compensa- 
tor order will be + 1. The only problem in determining the compensator 

order is dealing with cases in which the matrix N is (4.2.26) will be singu- 
lar. This will occur when some plant numerators contain identical zeros. 

This could be checked, but it essentially corresponds to a poorly posed prob- 
lem. 

4) Solution of Compensator Numerator Polynomials - 
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Once the plant numerator polynomials, the optimal loop gain polynomial, and 
the compensator order are found; a set of compensator poles are chosen and the 
compensator numerator polynomials are found by solving the mCn^+l) set of 

linear equations defined by (4.2.26). The solution of these equations depends 
on the condition of the matrix N, which is discussed in Section 4.4. Assuming 
that a solution does exist. We use an iterative solution routine. This is 
performed only once per design, and since the accuracy of the pole/zero can- 
cellations depends heavily on an accurate solution, the extra computational 
cost find storage required for the iterative solution is considered worthwhile. 


5) State-Space Realization of the LTR Compensator - 

Since the compensator is found in polynomial form it must be transformed to 
a state space form. I use a observable-canonical representation [CHI] which 
is as follows: 


0 0 -D- "I 



1 0 ,° 


1,0 m,0 

0 1 0 . 

, B = 

* * * * 

0 0 1 -D 


^ 4 « • » XI 



1,n N m,n N 


C = [ c 0 0] 

where is the ith coefficient of the compensator denominator polynomial, 
n^,j is the jth coefficient of the ith compensator numerator polynomial, and c 
is the product of the compensator poles which are not cancelled. 

Since the polynomial coefficients are sometimes very large, and will possi- 
bly contribute to later numerical problems, the state-space representation is 


134 



balanced so as to minimize the 1-norm of the matrix. This is a relatively 
cheap operation (computational cost proportional to n 2 ), and does result in 
better numerical properties for the robustness calculations. 
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